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The theory presented is based on a simple Hamiltonian for a vortex lattice in a weak impurity 
background which includes linear elasticity and plasticity, the latter in the form of integer valued 
fields accounting for defects. Within a quadratic approximation in the impurity potential, we find a 
first-order Bragg-glass, vortex-glass transition line showing a reentrant behavior for superconductors 
with a melting line near Hc2- Going beyond the quadratic approximation by using the variational 
approach of Mezard and Parisi established for random manifolds, we obtain a phase diagram con- 
taining a third-order glass transition line. The glass transition line separates the vortex glass and 
the vortex liquid. Furthermore, we find a unified first-order line consisting of the melting transition 
between the Bragg glass and the vortex liquid phase as well as a disorder induced first-order line 
between the Bragg glass and the vortex glass phase. The reentrant behavior of this line within the 
quadratic approach mentioned above vanished. We calculate the entropy and magnetic induction 
jumps over the first-order line. 

PACS numbers: 74.25.Qt, 74.72.Bk 



I. INTRODUCTION 

The phase diagram of high-Tj, superconductors in the 
H — T-plane is dominated by the interplay of thermal 
fluctuations and disorder [1, 2]. It is believed that at 
low magnetic fields near Tc the vortex solid melts into 
a vortex liquid (VL) via a first-order melting transition. 
Prominent examples of high-Tc superconductors exhibit- 
ing a solid-liquid melting are the anisotropic compound 
YBa2Cu307_5 (YBCO), and the strongly layered com- 
pound Bi2Sr2CaCu208 (BSCCO). When including weak 
pinning, the solid phase becomes a quasi-long-range or- 
dered Bragg glass (BG) [1] . At higher magnetic fields, the 
quasi-long-range order is destroyed and there exist also a 
vortex glass (VG) phase. The transition is marked by the 
disappearance of Bragg peaks in scattering data. There 
is strong experimental evidence especially for BSCCO 
[3, 4] but also for YBCO [5] that the BG-VG transition is 
first order, although in YBCO it has not been confirmed 
that this is really a proper phase transition, not just a 
crossover. So far, the transition line has been identified 
only by some magnetic anomalies in the response to the 
external magnetic field. 

For BSCCO it seems that the two melting lines arc 
part of a unified first-order transition line. For YBCO 
there are two possible experimental scenarios: First, the 
BG-VG and the BG-VL transition lines meet in a mul- 
ticritical point (MCP) [6] where the thermodynamical 
character of the BG-VG line is not clear. The first-order 
BG-VL melting line should continue beyond the multi- 
critical point ending in a critical point (CP) [7, 8] where 
a new fluid- like phase, the slush phase (VS), emerges (see 
Fig. 1). 

The second scenario consists of a unifled melting line 
[5] without a multicritical point which is the case for 
BSCCO. We note that the experimentally realized sce- 
nario is strongly sensitive on the doping of the supcrcon- 
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FIG. 1: Sketch of the two possible scenarios of the phase dia- 
gram of YBCO or similar high-Tc superconductors where the 
phase transition line lies near Hc2- The straight line corre- 
sponds to the BG-VG, BG-VL first-order lines with an exten- 
sion of the first-order line beyond the lower multicritical point 
(MCP). This is the first scenario discussed in the text. The 
dashed line corresponds to a unified BG-VG, BG-VL first- 
order line without a slush phase (VS) corresponding to the 
second scenario. The dashed-dotted line is the glass transi- 
tion line. The point GP is the intersection point of the glass 
transition line with the BG-VG, BG-VL line. 



ductor [8]. 

Beside the flrst-order transition lines, there exits a 
glass transition line between the VG and VS phases, if 
the VS phase exists, or between the VG and VL phases, 
if VS is absent. This glass transition line was predicted 
by Fisher et al. in Ref. 9, and observed when confirm- 
ing scaling rules for special current-voltage characteris- 
tics across the transition line [10]. Alternatively it was 
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proposed in Ref. 11 that the glass transition is window 
glass like with no scaling. Some people define an irre- 
versibility line beyond which magnetization sweeps are 
no longer reversible [12]. This seems to coincide with 
the glass transition line. A direct experimental determi- 
nation of the order of the glass transition in the vortex 
system of YBCO has not yet been possible. For BSCCO, 
there is recent experimental evidence that the glass tran- 
sition line could be of second order [13]. A sketch of the 
phase diagram for YBCO which contains the two scenar- 
ios is shown in Fig. 1. 

Beside experiments to determine the phase diagram, 
information comes from computer simulations based on 
the Langevin equation [14] or on frustrated XY-models 
[15-17]. The Langevin simulations confirm the second 
phase scenario without a slush phase, and the existence 
of a multicritical point on the melting line being unclear. 
In the frustrated XY model, the existence of a slush phase 
and of a multicritical point are also controversial [18]. In 
addition, Lidmar [19] carried out a Monte-Carlo simu- 
lation based on a defect model where he only obtains a 
first-order melting line and a glass transition line, but 
not a VS phase. 

Analytic approaches are based mainly on the Ginzburg 
Landau model [20], which is especially useful for YBCO, 
the cage model [21], or the elasticity model of the vor- 
tex lattice [22-28] with pinning. The Ginzburg Landau 
model with pinning was analyzed recently by Li et al. 
[20] where they found a phase diagram of the second 
scenario, with a single first-order melting line between 
the VG and VL phases as well as between the BG and 
VG phases, without an additional slush phase. The cal- 
culation was restricted to second order in the disorder 
potential. In a recent paper they also carry out an anal- 
ysis of a possible glass transition line in the fluid phase 
of the Ginzburg-Landau model where they found such 
a line only under a certain disorder model [29] by us- 
ing replica symmetry-breaking techniques which wc also 
use in this paper. In [21-23, 25-27] the phcnomcnologi- 
cal Lindemann criterion extended to include pinning was 
used in order to calculate the BG-VG and BG-VL transi- 
tion lines. In [24, 28] defects were taken into account for 
determining the transition lines. These approaches allow 
for an explanation of both phase scenarios. 

It is the purpose of this paper to investigate the above 
phase transitions in a defect melting model which was re- 
cently constructed for the study of defect-induced melt- 
ing of square (YBCO) and triangular (BSCCO) vortex 
lattices. The model is a modification of a simpler version 
in Ref. 30 which explained the melting transition of ordi- 
nary crystals by the statistical mechanics of defects on a 
hypothetical square lattice. This model was generalized 
for two-dimensional triangular crystals in Ref. 31. The 
model is Gaussian in the elastic strains and takes into 
account the defect degrees of freedom by integer valued 
gauge fields. The melting line is found from a lowest- 
order approximation, in which one identifies the melting 
point with the intersection of the high-temperature ex- 



pansion of the free energy density dominated by defect 
fluctuations with the low-temperature expansion domi- 
nated by elastic fluctuations. 

In this paper wc shall consider, in addition, the effect of 
weak disorder on the melting line near Hc2 ■ This will lead 
to a determination of the BG-VG and BG-VL transition 
line. The most prominent example for a high- Tc super- 
conductors with such a melting line is YBCO but also 
superconductors with a low critical temperature Tc such 
as BCS type superconductors or with a small anisotropy 
factor should have a melting line near Hc2- For concrete 
calculations, wc will restrict us in the following to the 
case of YBCO. 

The paper will first review the model and derive an 
effective Hamiltonian for the vortex lattice in the low- 
temperature solid phase and the high-temperature fluid 
phase without disorder. The model has two mutually rep- 
resentations. One can be evaluated efficiently in the low- 
temperature phase, the other in the high-temperature 
phase. The lowest approximation to the former contains 
only elastic fluctuations of the vortex lattice without de- 
fects. The dual representation sums over all integer- 
valued stress configurations, which to lowest approxima- 
tion arc completely frozen out. The tranverse part of 
the vortex fluctuations in the high-temperature approxi- 
mation corresponds to non-interacting three-dimensional 
elastic strings where the length in z-direction is dis- 
cretized with the dislocation length as the lattice spacing 
[32]. It is well known, that the lower critical dimension 
for an elastic string in a random potential [33] is three. 
This dimension separates the string system in higher di- 
mensions with two phases (a disorder dominated low- 
temperature phase and a temperature dominated high- 
temperature phase) from a single disorder dominated 
phase in lower dimensions. We encounter a similar situ- 
ation for the high temperature Hamiltonian in Sect. V. 
This is the reason, why we shall have to consider higher- 
order expansion terms [34]. 

We shall flrst expand the free energy to lowest order in 
the disorder potential in Section III. The result will be a 
unifled melting line. This line bends to lower magnetic 
fields in the direction for decreasing temperatures due to 
the disorder, in agreement with experiments. We obtain 
a remarkable reentrant behavior for this line. We do not 
obtain, however, a good agreement with experiments at 
low magnetic fields. In order to get better agreement with 
experiment and to determine also the glass transition 
line we further calculate, in the solid low-temperature 
phase and in the fluid high-temperature phase, the free 
energy non-pcrturbativcly by using once the replica-trick 
and further the variational approach set up by Mczard 
and Parisi [36] for random manifolds and spin-glasses 
[37]. It is based on replacing the non-quadratic part of 
the replicated Hamiltonian by quadratic one, with pos- 
sible mixing of replica fields. A transition line from a 
liquid to a glass consists within the Mezard-Parisi ap- 
proach on a boundary in thermodynamical space from 
a replica symmetric quadratic Hamiltonian to a Hamil- 
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tonian which breaks the symmetry in the rephca fields. 
The best quadratic Hamiltonian in the low-temperaturc 
solid phase is full replica symmetry broken correspond- 
ing to the BG-phase. In the high-temperature phase we 
find a region where the solution is full replica symmetric 
corresponding to the VL phase. Furthermore, we find 
a glassy region (VG) were the optimal quadratic Hamil- 
tonian depends on the form of the disorder correlation 
function. By carrying a comprehensive stability analysis 
in Section VII we show that for the form of the glassy 
state the kurtosis ki defined in (85) as functional on the 
positional disorder function is relevant. A Gaussian cor- 
relation function has kurtosis ki = 1. For high magnetic 
fields near Hc2 we obtain: 

For Ki < 1 we get a one-step replica symmetry broken 
solution with a third-order phase transition line. This 
corresponds to a correlation function with flatter tip and 
smaller tail than the Gaussian correlation function. In 
the case ki > 1 we obtain a full replica symmetry bro- 
ken solution. The free energy has the same form as in 
the one-step replica symmetry breaking case leading also 
to a third-order glass transition line. Disorder correla- 
tion functions with smaller tips and larger tails than the 
Gaussian correlation function belong to this case. For 
lower magnetic fields we obtain that the border in the 
disorder function space of one-step and continous replica 
symmetry broken solutions moves to lower kurtosis. 

The VG-VL phase transition happens just at the de- 
pinning temperature of a one-dimensional string in three 
dimensions subjected to impurities [1]. We calculate the 
free energies in the low-temperature solid and in both 
high-temperature phases. By the intersection criterion 
we obtain further the first-order BG-VL and BG-VG 
transition line. We do not find an additional first-order 
line which separates the slush phase VS from the vortex 
liquid VL. Summarizing, we obtain within our approach 
only the second scenario. In the low-temperature solid 
phase our analysis corresponds to the analysis of Kor- 
shunov [38] and Giamarchi and Doussal [39] using the 
Mezard-Parisi approach for the vortex lattice system in 
random potentials. Because this system does not con- 
tain defects only the Bragg glass phase can be described 
correctly. Within our approach we include beside the dis- 
order ones also the defect degrees of freedom by integer 
valued fields which arc important to obtain the melting 
transition. This allows us to compute the whole phase 
diagram for YBCO. In this sense our theory is a direct 
generalization of the earlier vortex lattice approaches in 
random potentials. 

Beside the results above, we give in Appendix B a 



derivation of some stability theorems of saddle point so- 
lutions similar to the theorems of Carlucci et al. [40] de- 
rived within the large A^'-limit approach of Mezard and 
Parisi [36] where N' are the number of components of the 
random manifold. Because in general the number of com- 
ponents of a given random manifold is rather small it is 
useful to generalize these theorems also to the variational 
approach of Mezard-Parisi not existent in literature yet. 

The paper is organized as follows: 
In Section II we state the model of the vortex lattice with 
defects and impurity degrees of freedoms. We derive in 
Section III the effective low and high-temperature Hamil- 
tonian of the vortex lattice without impurities. With the 
help of these effective Hamiltonians we calculate in Sec- 
tion IV the BG-VG, BG-VL transition line within the 
second order perturbation theory in the disorder poten- 
tial. In Section V we introduce the Mezard-Parisi varia- 
tional approach. Section VI calculates the saddle point 
solutions of the self-energy matrices for the variational 
free energy within the Mezard-Parisi approach in the 
fiuid high-temperature phase. Section VII deals with the 
stability of the calculated saddle-point solutions. Sec- 
tion VIII calculates the saddle point solutions in the solid 
phase. In Section IX, we discuss the phase diagram of the 
Mezard-Parisi approach for YBCO and compare it with 
the experimental ones. Furthermore jump quantities are 
calculated in this section. Section X contains a summary 
of the paper. 



II. MODEL 

The partition function used here for the vortex lattice 
without disorder was proposed in Rcf. 32. It is motivated 
by similar melting models for two-dimensional square [30] 
and triangular [31] crystals. Motivated by the fact that 
YBCO has a square vortex lattice we restrict us here to 
a discussion of the phase diagram of such type of lat- 
tice. The generalization to triangular vortex lattices is 
straight forward [32] resulting only in a slight difference 
in numerical values. We briefly summarize the important 
features of the model. The partition function of the dis- 
ordered flux line lattice can be written in the canonical 
form as a functional integral 

Zfl^l 2?[^.„a.™,n,]e-(^"["--™-"'l+^-["-l)/'^--^, (1) 

where 
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El ^ 2 , 1 2 f'^'^i \ Cll — 2C66 /sr^^i \ , '^66 



27rz ^ [ 



(2) 



is the canonical representation of elastic and plastic 
energies summed over the lattice sites x of a three- 
dimensional lattice, and where (T21 = o'i2 are stress 
fields which are canonically conjugate to the distortion 
fields [30]. The subscripts i,j have the values 1,2, and 
l,m,n run from 1 to 3. The parameter /3 is proportional 
to the inverse temperature, /3 = U'^cqq / k bT {2t:Y , 
where a is the transverse distance of neighboring vortex 
lines, and is the persistence length along the disloca- 
tion lines introduced in Ref. 32. Note that 03 is assumed 
to be independent on the disorder potential on the av- 
erage. The volume of the fundamental cell v is equal to 
a^oa for the square lattice. 

The matrix Nij{x.) in Eq. (3) is a discrete- valued local 
defect matrix composed of integer-valued defect gauge 
fields 7ii,7i2. It depends on the lattice symmetry [31]. 
For a square vortex lattice it is given by 



ni 712 
7i2 -ni 



(3) 



The lattice derivatives Vm and their conjugate coun- 
terparts Vm are the lattice differences for a cubic three- 
dimensional crystal. In the xy-plane they are defined by 



VJ(x) = [/(x)-/(x-ae,)]/a, 
VJ(x) = [/(x + ae,)-/(x)]/a 



(4) 



for a lattice function /(x), where are unit vectors 
to the nearest neighbors in the plane. The correspond- 
ing derivatives in z-direction are defined similarly. We 
have suppressed the spatial arguments of the elasticity 
parameters, which are functional matrices Cy (x, x') = 
Cij{x. — x'). Their precise forms were first calculated by 
Brandt [41] and generalized in Ref. [32] by taking into 
account thermal softening relevant for BSCCO. 
The second term in the exponent of (1) 



Hdis[ui] = ^y(x-l-u). 



(5) 



accounts for disorder. The measure of the functional in- 
tegral is 



I?[u,;,(T,>„,ni]=det 



C66 


1/2 

det 


1 


4(cii-C66). 


_27r/3_ 





n 



/oo 
da. 



n E 

. j nj(x) = — oc-J L 



-| 5/2 

du 

a 



(6) 



The disorder potential V{x.) due to pinning is assumed 
to possess the Gaussian short-scale correlation function 



T/(x)y(x') = A{xi~x'i)d, 



(7) 



d(T) 03 



where K{xi — x-) w 1/(C')^ for |x — x'| < and is 
zero elsewhere, and 0o is the magnetic flux quantum 
00 ~ hc/2e. The parameter ^' is the correlation length 
of the impurity potential which is similar to the coher- 
ence length S^ab in the xj/-plane. \ab = Aafc/(1 — b) is the 
screened penetration depth in the xy-plane [23] . 

The temperature dependence of the parameter d{T) is 
mainly due to the temperature dependence of the corre- 
lation length and the pinning mechanism where we dis- 
cuss in the following the JTc-pinning or (5;-pinning mech- 
anisms [1]. 

Both pinning mechanisms are extensively discussed in 
the review of Blatter et al. [1]. We just mention that the 
(5Tc-pinning mechanism has its origin in fluctuations in 
in the Ginzburg-Landau free energy and the (5Z-pinning 
mechanism is due to fluctuations in the mean free path 
coming from fluctuations in the impurity density. The 
parameter d{T) is different for both pinning mechanisms 
[1]: 

d{T) = do(l - T/Tc)-^/^ for 6T, - pinning , (8) 
d{T) = do{l-T/Tcf/^ for ,5? - pinning. (9) 

The correlation functions for both mechanisms can be 
derived in Fourier space by taking into account the or- 
der parameter shape of a single vortex [1]. This is of 
long range, resulting in a divergence of the Fourier trans- 
formed disorder correlation function K{q) at q = for 
the 5Tc-pinning mechanism. This divergence is regu- 
lated for a vortex in a lattice by omitting the regime 
g < 1/a because the order parameter of the superposi- 
tion of non-cut off single vortex order parameters on the 
lattice would otherwise scale with the system size. We 
shall see below in Section VII that this is the momentum 
region of the disorder correlation function which deter- 
mines mainly the form of the free energy in the fluid 
phase near the glass transition line and thus the order 
of the glass transition. Other correlation mechanism as 
for example screening of impurities are not taken into 
account in these single vortex disorder correlation func- 
tions. The screening of impurity potentials is important 
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because the nearest neighbor distance between impuri- 
ties is typicaUy of the same size as the coherence length 
U [!]• 

AU this leads us to use in the calculations to follow an 
effective disorder correlation function with the Fourier 
transform 



K{q) = 27rexp(- 



/2) 



(10) 



leading also to an exponentially vanishing of the disor- 
der correlation function in real space. The advantage for 
using this effective correlation function is that one gets 
simple analytical formulas in the calculations. The pa- 
rameter ^' in (10) is an effective correlation length which 
can also include for example screening effects of the im- 
purities in the 5i-pinning case. The approximation (10) 
leads to well known approximations for the temperature 
softening of quantities which use the disorder correlation 
functions as an input as for example the temperature 
softening of the coherently time-averaged pinning ener- 
gies [1]. 

In the following sections, we come back to the more 
general case without the assumption (10) for the disor- 
der correlation function especially in Section VII where 
we show that the order of the glass transition line depends 
strongly on the form of the correlation function. The free 
high-temperature energy formulas (16), (55), (64), (75) 
and (100); are valid irrespective of the form of the disor- 
der correlation function. In the low-temperature regime, 
the form of the energy expression cannot be find out for 
general correlation functions. In this case we restrict us 
to the effective disorder potential (10) where in contrast 
to the glass transition line, the final energy expression 
should not change much when changing the disorder po- 
tential. 



III. PARTITION FUNCTION OF SOLID AND 
FLUID PHASES FOR V = 



In this Section we determine the partition function of 
the low temperature phase (solid phase) and the high- 
temperature phase (fluid phase) for ^ ~ 0. This was 
done before in Ref. 32. Here, we give similar expressions 
which are appropriate for calculating correlation func- 
tions of vortex displacements useful for a discussion of 
the disorder problem. 

For the low-temperature limit of the partition function 
in (3) we first integrate out the stress fields (Jij. Then 
the low-temperature part of the partition function corre- 
sponds to the defect configuration rii = 0. This results 
in a fluctuating part of the form 



\ r «.(x)i 




7-00 a . 


cxp 





1 



with the low-temperature Hamiltonian 

Hf)[Ui] ^ Ht^oIUi] = I ^(VjMj)(cil - 2c66)(ViUi) 



-(ViWj -I- VjUi) C66 (ViWj -I- VjUi) + (Vsui) C44 (Vaui) 



= ^ ^i^iUL) cii (Viiii) + (VsUl) C44 (V3UL) 

X 

+ {ViUr) C66 (ViUr) + (Vswt) C44 (Vsut) 



(12) 



and the normalization factor J\f = 1. Here ul = Plu 
is the longitudinal part of the displacements where the 
proje ctor Pl is given by {PL)jk ^ -(Wl^fD^j <^ 
(l/y|V|I)Vfc. The transversal part of the displacements 
is then given by ut = Ptu = u — u^. The corrections 
to the fluctuating part of the free energy — ln(Zfl)/fcBT 
in the low-temperature expansion is exponentially van- 
ishing with an exponent proportional to —l/ksT [30]. 

For the high-temperature limit of the partition func- 
tion (1) we carry out first the sum over the defect fields 
rii, n2. By a redefinition of the stress fields (jg = 
(ail + (T22) and (T„ = (an — (T22) we obtain that ai2 
and au can only have integer numbers. The lowest-order 
terms in the high-temperature expansion of the partition 
function (1) for V = corresponds to au = CTm = 0. 
After carrying out the integrals over the stress fields Ug 
and (Ti3 we obtain a partition function of the form (11) 
with 

Ho[Ui\ = i?T-»ooK] = ^'^{^2Ui){cii - C66)(ViUi) 

X 

+ (Vaui) C44 (Vaui) 

= 7;'^(^i'^L) (Cll - C66)(VjUl) + (VsMi) C44 (VsUl) 



+ (VaitT) C44 (VsitT) 



(13) 



(11) 



and Af ~ l/(47r/3)^. Similar as in the case of the low- 
temperature expansion one can show that corrections to 
the fluctuating part of the free energy — ln(Zfl)//cBr due 
to non-zero terms in the stress fields ai2 ^ or (t„ ^ 
in the high-temperature expansion are exponentially 
vanishing with an exponent proportional to —ksT [30]. 

Expression (13) shows the remarkable fact that the 
transverse part of the high-temperature Hamiltonian (13) 
is effectively one-dimensional with a non-zero dispersion 
only in z-direction. This results in diverging thermal fluc- 
tuations in Ut- In contrast to this we obtain for that 
part of the Hamiltonian (13) corresponding to longitudi- 
nal fluctuations an effectively three-dimensional Hamil- 
tonian as in the low-temperature case (12) with finite 
temperature fluctuations. This can be better understood 
by the fact that only the transverse fluctuating part of 
the vortices couples to the defect fields while the longitu- 
dinal part is still not effected by them [32, 42, 43]. The 
reason is that the flux lines in a vortex lattice cannot 
be broken which means that defect lines are confined in 
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the plane spanned by their Burger's vector and the mag- 
netic field. In conventional crystals we do not have such 
a constrained [30]. It then clear that the large thermal 
fluctuations of the transverse part results in a destruction 
of the long range order in the sense that Bragg peaks are 
vanishing in the fluid phase. 

Summarizing, with the help of the stress representation 
(2) we obtained the lowest-order Hamiltonians for the 
solid (12) and the fluid phase (13). We saw further that 
the higher-order corrections to this lowest-order results 
corresponds to integer-valued defect contributions ^ 
in the solid phase, signals for the liquid, and integer 
valued stress contributions (7i2 7^ or (t„ in the 
fluid phase, signals typical for a solid. In the following 
we restrict us to the lowest-order Hamiltonians (12) for 
the solid phase and (13) for the fluid phase to discuss 
disorder corrections in both phases. 



disorder constants 
Do{2{u^)) = d{T) 



as 



D^{q)=d{T) 



«3 d^Ul^K^Q) / 



{ksTY 2^ / d(g2/2) 



k{q) e-^<"'>, 
(17) 

(18) 



where V{2{u'^)) = Dq{2{u^)). Note that we have X'(O) = 
2?oo(0) = I?o(0) for the Gaussian correlation function 
(10). For this correlation function, we obtain 



V{2{u')) « d{T) 



as 00 ^ab 



{ksTy 



■^ab 



(19) 



Furthermore, we define the corresponding disorder cor- 
relation lengths by 



IV. QUADRATIC APPROXIMATION IN 
DISORDER STRENGTH 



To lowest non- vanishing order in the disorder potential 
V we obtain for the first non-vanishing term in the free 
energy F ~ —kBThi{Z) a term proportional T^^ given 
by 



1 



fl.V2 



2kBT 



(14) 



J2 (^(x + u)F(x' + u)) - (T/(x + u)) (y (x' + u)) 



Note that the dimension of A is {ksT)'^ (7). 

We restrict us to the diagonal summands x = x' 
where non-diagonal terms results in corrections only to 
the low-temperature expansion of -Ffl.v'2 being a factor 
((m^) + ^'^y/^/a smaller than the diagonal terms. The 
calculation can be most easily done by working in the 
Fourier representation 

By using (12) and (13) we obtain for the low- and high- 
temperature part of the free energy 



-N 



d^ 



q 



2iksT) J (2.)3-^(^)^(-^) 
X (1 - cxp [-q- (uu)t^o ■ q]) , 

— N f d^n 

^/(4^(^)^(-^)- 



(15) 
(16) 



For the determination of the melting line only the sec- 
ond term in the bracket of F^y2 is relevant because of 
a cancellation when determining the intersection of the 
high and low-temperature expansions of the free energy. 

After carrying out the momentum integral and disor- 
der averaging we obtain for this term N'DkBT/2 with the 
disorder constant P defined by the help of the generalized 
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772 
s 

1 



1 



(27r) 
A-(0) 



d\K{q) im. 



d 



d{q^li^ 



■K{0) 



(20) 
(21) 



S 00 

with ^'^ = = Co for the Gaussian correlation func- 
tion (10). In the following, we carry out the calculation 
of the free energies explicitly for the Gaussian correla- 
tion function where we use V and ^' without indices. As 
mentioned in the introduction our final results in this 
section and also for the fiuid phase in the Mezard-Parisi 
approach are more general valid without restrictions on 
the disorder correlation functions. 

By recalling the results for Zfl without disorder [32] 
with the low-temperature Hamiltonian (12) we obtain 



2N 



1 



det [(27r/3)c44/c66] 

and for the high-temperature part (13) 
'a3\2W 1 1 



r^l -too 



£{1,6} 



hi 



a J 2^ det [(27r/3)2c44/c66] 



(22) 



(23) 



with 



1 1 

2Vb^ 



2Vbz 



d kdk^ In 



BZ 



^K*K, 
C44 ^ 



d kdk-i In 



BZ 



1 



Cll 



C66 1^*3^3 



C44 



(24) 



where K,n is the eigenvalue of iVm. The k^k^- 
integrations in (22) run over the BriouUin zone of the 
vortex lattice of volume Fez = {2t:Y /v. According to 
the intersection criterion we equate (22) and (23) and 
obtain the equation for the temperature [32] 



kBT 



1 



V det 



l/N 



C66\ 



(25) 
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FIG. 2: Unified BG-VG, VG-VL transition line Bm{T) (29) as 
a function of the temperature. The curves in the upper figure 
are calculated for ^Tc-pinning mechanism (8), the lower for 5;- 
pinning mechanism (9). The parameter do of the solid curves 
are chosen such that we get a good fitting to the experimen- 
tally determined VG-VL line by Bouquet et al. [6] (dashed 
curve with square points) for both pinning mechanisms. The 
dotted curves are variations from these best fitting curves 
given by disorder parameters = (1 ± 1/2) do where do are 
the disorder parameters of the solid curves of both mecha- 
nisms given by 27rdo^a6/^' ~ 8.5 • 10~* {STc- pinning) and 
27rdo Cab/^'^ = 1-01 ■ 10"*^ {51- pinning). The solid curves with 
the stars are calculated by solving (25) with elastic moduli in 
the range 6 < 0.2 with do given above. 



The solution determines the first-order BG-VG, BG- 
VL transition line with disorder. The disorder enters 
the equation via the disorder function T). Analytic ex- 
pressions can be obtained by taking into account that 
C66)C44 cii. This implies that we can neglect h and 
^11 in (22). 

Brandt [41] determined the elastic constants for two 
different regimes b < 0.2 and h > 0.5 where b = 
B/Hc2{T). We shall see below that for YBCO we have to 
determine (25) in both regimes to find the entire relevant 
part of the BG-VG and BG-VL line. The most important 
part, however, lies in the regime b > 0.5 which will now 
be treated explicitly. In this regime the elastic moduli 
C44 and C66 are given by [41] 



C66 = 0.71 (1 - b) 



C44 



47r(l + A2fc2 + A2,fc|) 167r2A2 



(26) 
(27) 



A = A/(l — &)"'^/^ is the screened penetration depth calcu- 
lated from the penetration depth A. In xy-dircction, we 
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FIG. 3: Logarithmic plot of the BG-VG, BG-VL first-order 
line Bm{T) (29). The experimental points (triangles) corre- 
spond to the experiment of Pal et al. [7] showing a reen- 
trant BG-VG line. The theoretical determined solid curves 
are derived from (29). The pure solid curve {b > 0.5) and the 
solid curve with the stars {b < 0.2) correspond to the lower 
curves in Fig. 2 (5Z-pinning). The theoretical curves are cal- 
culated with the disorder parameter 27rdo Cat/?' = 1-01 ' 10~® 



of Fig. 2. 



denote it by Xab, and in z-direction by Ac. For YBCO 
we have [35] A(T) = A(0)[1 - {T/T,)]-^/^ UiT) = 

umi - (T/r,)]-V2. 

For later use, we define the Lindcmann parameter [32] 



a^v Vbz Jbz 



f d^kdks— V : 

Jbz C44 211^ 



1 



i=l,6 " 



kaT, 



C44(^,0)C66(^,0) 



1/2 



(28) 



given by c\ = 



>o/a^ where the average is taken 



with respect to the low-temperature Hamiltonian (12) 
representing the elastic energy of the vortex lattice. K^z 
is the boundary of the circular Brillouin zone K^^, = 
47rB/0o. For YBCO, we obtain [32] cl « 0.18 on the 
melting line without disorder in accordance with typi- 
cal Lindcmann numbers for crystals [30]. We note that 
this number does not depend on the magnetic field which 
specifies the point on the melting line. In the following 
we denote C44 (Xbz/\/2,0) and C66(i^Bz/%/2, 0) in final 
expressions as for example in (28) by the abbreviations 
C44 and C66. 

From (25) we can easily calculate the BG-VG, VG-VL 
line. By taking into account c^a^ ^ ^'^ which results 
in V{2{u'^)) « X'(O) for YBCO we obtain for the unified 
BG-VG, VG-VL line 



0^(1-6)3 3.9-10- 



(29) 



Here we used a — \f^l\lB and the typical defect length 
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as ~ 4a%/2 Aat/Ac\/7r(l 
disorder function 



bf/'^ [32], which results in the 



2?o(0) 



So 



{kBTf 



(30) 



Note that (29) is valid irrespective of the disorder corre- 
lation function. 

Parameter values for optimally doped YBCO were 
given by Ref. 35 as Aah(O) w 1186A, ^^6(0) ~ 15A, 
Tc ~ 92.7K. The Cu02 double layer spacing is Us — 12A, 
and the anisotropy parameter 7 = Xc/Xab is approxi- 
mately equal to 5. From (30) we obtain a unified BG-VG, 
BG-VL hue which scales like ~ e"^^ /^^. Here Ay 
is some constant independent of Bm- This results in a 
BG-VG, BG-VL line showing a reentrant behavior. In 
Fig. 2 we show B^n with i5Tc-pinning on the upper figure 
and (5Z-pinning on the lower figure for various values do- 
do of the straight line curves is chosen such that we have 
approximately the best experimental curve fitting to the 
BG-VG, BG-VL curve of Bouquet et al. [6] shown by 
the dashed line with square points. We obtain in fact a 
reentrant behavior of the BG-VG, BG-VL line. This is in 
accordance with the quadratic in disorder calculation for 
YBCO within the Ginzburg-Landau approach by Li et al. 
[20]. A reentrant behavior of the BG-VG line was also 
seen in the experiments of Pal et al. [7] and Stamopoulos 
et al. [44]. We must clarify that these experiments are 
in contradiction to the majority of experiments which do 
not see any reentrant behavior of the BG-VG line (see for 
example [5, 6, 12, 45, 46]). The discrepancy in the shape 
of the BG-VG transition line lies presumably in differ- 
ent physical setting of these experiments to the standard 
ones showing no reentrant behavior. The experiment of 
Pal et al. uses a crystal with a low density of twins which 
could lead to deviations in the shape of the BG-VG line 
due to [47] . The experiments of Stamopoulos et al. mea- 
sures the ac permeabilities which drives the crystal out 
of thermodynamical equilibrium. 

The solid lines with the stars at the left hand side of 
Fig. 2 are calculated by solving (25) restricted to the 
transverse fluctuations with the elastic moduli in the 
range b < 0.2 given in [32, 41]. To calculate the tran- 
sition curves with the elastic moduli b < 0.2 as well is 
rather important because the solid curves in Fig. 2 cal- 
culated with moduli b > 0.5 reaches immediately the 
range b < 0.2. Note that the solid curves and the solid 
curves with stars are calculated by using the same disor- 
der constant do of values 27rdo ^^b/^'^ = 8-5 • 10"^ {STc- 
pinning) and 27rdo ^^f,/^'^ = 1-01 ' 10"^ ((5/-pinning). 

We obtain from Fig. 2 that the curves of the Sl-pmning 
mechanism fits much better to the experimentally given 
BG-VG, BG-VL hue than the 5Tc-pinning curves. This 
is in accordance to the observation in Ref. 46. Whatever 
the correct experimental BG-VG, BG-VL line shows a 
reentrant behavior or not, it is not satisfactory within 
our approach which is restricted to second order in the 
impurity potential, that the solid curves with the stars 



(6 < 0.2) goes to zero at T Ri Tc. That this is true can 
be best seen in a logarithmic plot of the BG-VG, BG-VL 
transition line B^ (29) which is shown in Fig. 3. The 
straight line and the straight line with the stars corre- 
spond to the curves of the (5Z-pinning mechanism shown 
in the lower part in Fig. 2. The dashed curve with the 
triangle points is the BG-VG, BG-VL line measured by 
Pal et al. in Ref. 7 mentioned above. This curve shows a 
reentrant behavior. Both curves are in disagreement at 
small magnetic fields. Thus, we should go beyond second 
order in the disorder strength to get better accordance 
with the experiments. This will be done in the following 
sections. 



V. REPLICA VARIATIONAL METHOD OF 
MEZARD AND PARISI 

In order to go beyond second-order perturbation the- 
ory in the impurity potential, we use the well known 
replica trick [48] InZ = lim„^o(l/")('^" ^ 1) where the 
overline means disorder averaging and 



X e 



-S„(ifo[<,o-?„,n^H-ffdiB[<])/fcBT 



X e 



(31) 



with 



x,x' 

(32) 

Here, the extra ksT term in the denominator in (32) 
comes from the disorder average. 

The average Z" cannot be calculated without further 
approximations. In the following we use the low- and 
high-temperature approximations of Section HI for the 
result after the integration over the stress fields ct," and 
defect fields nf in (31). Thus, we have to calculate 
partition functions of the following form 



z« = n^n 

a x,i 

X e-^/^-^, 
with the total Hamiltonian 



<(x) 



(33) 



(34) 



where Hq is given by (12) in the solid phase and by 
(13) in the fluid phase. Both are complicated expres- 
sions which will need further approximations. The com- 
plications comes from the large replica mixing interaction 
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part iJjjg in (31). In the following, we shall use a vari- 
ational replica method which was first given by Mezard 
and Parisi [36] . With the help of this method also used 
before for random spin models [37] they were able to cal- 
culate the glass transitions of isotropic random manifold 
systems. These systems are described by the Hamilto- 



Hrm = / d^-'^x 



[-u(x)(V. V)u(x) + nu)]. (35) 



Here u is an A'^'-dimensional vector describing an 
A'^'-dimensional manifold embedded in a d-dimcnsional 
space. V is an impurity potential with a certain corre- 
lation function. When comparing the solid Hamiltonian 
(12) with the random manifold Hamiltonian (35) and fur- 
ther by setting the correlation length ^' = in (19) we 
obtain that the transversal part of (12) looks similar to 
a random manifold with d — iV' = 3in(i = 5 dimensions 
[36] with a delta-like impurity correlation potential. In 
the fluid phase described by the high-temperature Hamil- 
tonian (13) we obtain for the transversal part a random 
manifold with d — N' = 1 and c? = 3 well known as a 
string embedded in three dimensions. The difference to 
the random manifold system comes then mainly from the 
discretization in the third direction by the dislocation 
length 03 relevant in the fluid phase as will be shown 
below. It is well known that there exist for N' > 2 
in a d — N' + 1 random manifold system correspond- 
ing to a string in d dimensions a roughening transition 
separating a low-temperature disorder dominated phase 
from a high-temperature thermal phase [33] . For N' < 2 
this phase transition is not existent and the system is 
dominated mainly by disorder fluctuations. It is now 
believed by computer simulations that at the critical di- 
mension N' = 2 corresponding to a string in three dimen- 
sions with a (5-correlated impurity potential the rough- 
ening transition of the string system is described by a 
crossover [33] . Below we show that this roughening tran- 
sition corresponds to the glass transition of the vortex 
lattice. That the vortex lattice at d = 3 is in fact at the 
lower critical dimension for a glass transition was men- 
tioned before for an XY-model of the gauge class type 
[49]. This XY-model as similar ones with other disorder 
potentials mentioned in the introduction are toy model 
for a disordered vortex lattice in superconductors. 

The Mezard-Parisi theory consists in replacing the 
non-quadratic part of the Hamiltonian as quadratic with 
a possible mixing of replica fields. By using the Bogoli- 
ubov variational principle we can find the best matrix of 
this quadratic form so that the free energy of the vari- 
ational Hamiltonian is as close as possible to the actual 
free energy of the system. This means that we have to 
search the minimum of the variational free energy 



with the harmonic trial Hamiltonian 

Atrial = I E E WG-^(X - X')U^(X') . 
x,x' a,/3 



(36) 



(37) 



Here (•)triai stands for the averaging with respect of the 
Gibb's measure of the trial Hamiltonian iJtriaii ^md -Ftriai 
denotes the associated free energy. In Section IX, we 
shall use the intersection criterion with i^var for the solid 
and fiuid phase to determine the BG-VG, BG-VL tran- 
sition line. By using (37), we obtain the free energy as- 
sociated with (33): 



-kuT- 



N 



d''kdk-i In 



2Vbz Jbz 
Tr{[G-i(k)-Go-i(k)l]G(k)} 

kBT\n{N) 



det f^G 




(38) 



trial 



where we use bold symbols for vectors and matrices in 
the vortex displacement plane. The symbol I denotes the 
unit matrix in replica space. The trace Tr[ . . . ] runs over 
the replica indices and vortex displacement indices. In 
principle we can obtain a general expression for the dis- 
order term given by the last term in (38). Because one 
should use different approximations for the solid phase 
and the liquid phase, we shall give directly approxima- 
tions for this term in both phases at the beginning of the 
following sections. 

It will be clear soon for the solid as well as the fluid 
phase that Gap can be chosen to have the form 



Q/3 



(39) 



where I is the two-dimensional unit matrix in the vor- 
tex displacement plane. To flnd a local minimum of 
(38) in the space of all symmetric self-energy matrices 
aap was simplified considerably by Parisi in the case of 
spin glasses. There he restricted the search of the mini- 
mum for (36) to the case of some sort of closed algebra 
known as the algebra of Parisi matrices [37, 50] . In Ap- 
pendix B we prove some stability theorems for stationary 
points of i^var (38). These are summarized at the end of 
Appendix B2. The restriction of the minimum search 
to self-energy matrices in the Parisi-algcbra is justified 
among others by the fact that a local minimum within 
the Parisi-algcbra is automatically a local minimum in 
the whole self-energy space without the restriction to the 
Parisi- algebra. This is shown in Appendix B. 

In general the minimum self-energy matrix cTq,^ is not 
symmetric under the interchange of replica indices which 
means that the local minimum aap of -FVar (38) is not 
unique. This is typical for glasses where the minimum of 
the free energy is degenerate [37]. This degeneracy cor- 
responds to the degeneracy of the stable states in glasses 
with high energy barriers between them. These are re- 
sponsible for the irreversibility phenomena beyond the 
glass transition lines in high-temperature superconduc- 
tors mentioned in the introduction. 
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VI. FLUID PHASE 



Gaa- This results in 



In this section, we derive the variational free energy 
i^var in the liquid phase. We obtain 



N 




trial 



d?q A{q)e 

a,l3 



G„„(0)+G3^j(0)-G„,5(0)-G3„(0) 



(40) 



with 



= ^tr [G„a(0) + G/50(O) - G„^(0) - G^„(0)] . 

(41) 

where A(q) is the two-dimensional Fourier transform of 
A(x). The trace tr[. . .] runs over the vortex displace- 
ment indices. In (40) we restricted us in the double sum 
over X, x' on the diagonal summands x = x'. The reason 
for the validity of this restriction comes from the obser- 
vation that due to (13) the non-diagonal summands are 
given by 



trial. 



(42) 



X e"5q-[GL(0)+Gj^,(0)-G;[,j(x-x')-G|-„(x-x')]q 
X e-5q-[GL(0)+G^,5(0)-G^^(x-x')-G§„(x-x')]q 

where G^ = PlGPl and G^ = PtGPt are the lon- 
gitudinal and transversal components of the Green func- 
tion. By using (13) the second exponent in (42) corre- 
sponding to the transversal part of the vortex fluctua- 
tions can be tranformed to 

1 

^ 2 



(q'')'(o) [GlM + Glpm 

(q^)2(x-x') [GI^(0) + G^J0)] 



(43) 



where G^^ (x) = Tr[G^^ (x)]. Due to the large thermal 
effective one-dimensional transverse fluctuations we have 
either G'^^(O) oo where Gl^[Q) - G^^(O) is flnite and 
a, j3 is arbitrary or G^^(O) oo and G^^(O) is flnite 
for a ^ (3 where in both cases the self-energy matrix is 
restricted to the Parisi algebra. This is shown in Ref. 36. 
From this we obtain the vanishing of (42). 

First, we take the variation of the free energy (38) with 
respect to the diagonal Green function matrix elements 



Uaf) = 

/3 



(44) 



That the minimum of Fvar should be found in the sym- 
metric self-energy matrices with the constraint (44) is 
suggestive because (44) justifles that the Hamiltonian 
(37) has the global translational symmetry u"(x) 
u"(x) -I- 1 for any vector which has also the disorder 
Hamiltonian (32). 

In the most general case within the Parisi-algebra, the 
form of the self-energy (Jap with the constraint (44) can 
be described by a continuous function (t{s) with < s < 
1 [36] . In that case the trial free energy takes the form 



A/ 



,ar - 4 lim - [F,„(B[A]) - F,ar(0)] 



N n^o n 







f ds 




'o 





fA(s) 



/var(O) = ^lim„^o-^'var(0) = -fcfiT 

iV n 



\{t^ I d^kdh In 



det 



2TTkBT^ 



(45) 



(46) 
Vo{0) 



where 
.9(A) 



1 



d kdk^Ti 



BZ 



1 Gn 



AI 



(47) 



The gap function A(s) and the self-energy function (t(s) 
corresponding to the self-energy matrix aap in the non- 
continuous case is related by 



, , dcr{s') 
^^^^ 



(48) 



B[A{m)] corresponding to Bap (41) in the continuous 
case is given by 



B[A(.)] = M:i,[A(,)] 



(49) 

In order to find the local minimum of /var we have to 
take the derivative of (45) with respect to A(m). This 
results in 



kuT 



V',{2B[A{s)]) , 



(50) 



where D'{x) is the derivative {d/ dx)D{x). We point out 
that (50) shows that 



cr(s) > 



A(s) > 0. 



(51) 



In the following, we discuss solutions of this equation in 
the case that (j{s) does not break the replica symmetry, is 
one-step replica symmetry-breaking or continuous replica 
symmetry-breaking . 
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A. Symmetric solution 

We now solve Eq. (50) for cr(s) with an Ansatz which 
does not break any replica symmetry. The Ansatz for 
a{s) in this case is 

<t{s) = ao . (52) 
By using the equations (48) and (50) we obtain 



a(s) = 0, 
A(s) = 0. 



(53) 
(54) 



From this we obtain that B[A{m)] oo for infinite area 
of the system. This results in 



A/v, 



0. 



(55) 



B. One-Step replica symmetry-breaking 



The simplest possible extension of the replica symmet- 
ric case above consists of a one-step replica symmetric 
solution given by 



fJo for < s < nil , 
(Ti for mi < s < 1 . 



(56) 



By using this Ansatz in (45) we obtain 




/ A, 

4 arcsinh — - 



2 -(l-mi)P(2^.g(Ai: 



(57) 



where we used = which can be derived from (50) and 
(48) similar to the replica symmetric case. Furthermore, 
we used the abbreviation Ai = nixipx — <7q) = miui 
and Ki = Aia§/c44. We restricted us in the calculation 
of A/var to the transversal components of Go which ef- 
fectively means that we set Gq = in the calculation 
g(A) in (47). The longitudinal term in .g(A) is a factor 
a^C44/aicii = (a^C44/a|c66)(c66/cii) ~ (7r/4) (cee/cn) 
smaller which is justified by cge ^ cn [41] irrespective of 
the value of Ai. For the derivation of A/var we used for 
g{A) in (47) 



3(A) 



1 



C44 VbZ 
1 



BZ [2 - 

1 a? 



2 cos(A:3a3) 



2 Al/2(1 + A/4)1/2C44 



(58) 



We now determine the stationary point of A/var (57) . By 
setting the derivative of A/var with respect to Ai and mi 



equal to zero we obtain two equations for the stationary 
values of Ai and mi. These are given by 



1 



1 



1 

mi 



.-1/2 



1 



(1 - mi) — 



^3/2 C44a2 



(59) 



1 1 



/ A, 

4arcsinh — - 



Ai 



1/2- 



H-Ai/4 



(60) 



In the following solution of (59) and (60) we use that 
Ai ^ 1 in the interesting range near the glass transition 
line which we expect at Ai = 0. This will be shown 
below. These two equations can be solved exactly in this 
limit resulting in 



ni\ = (r'(o)A)- 



a; 



1/2 



2A~ 



- 1 



mi 



(61) 
(62) 



where constant A similar to the Lindemann constant 
written for general disorder correlation functions (see the 
definitions (17)-(21)) 



^0,00 = 



4 C44a^ 



0,00 



0.3 



(63) 



with A = Aq = Aoo for the Gaussian correlation function 
(10). Here, we mention that A « 6/27rc| ^ 1 near the 
melting line without disorder V — [32]. This is the 
magnetic-temperature regime, we are interested in. By 
the help of (57), (61) and (62) we can calculate the free 
energy A/var getting 



Af ^M^p 



(0) [l-(I?^(0)A^)-i/3 



(64) 



for X'oo(0)Aoo > 1 in the regime {{V{0)A)^^^ - 1)/A < 1 
for the Gaussian correlation function (10). As suggested 
by the indices, expression (64) is more general valid ir- 
respective of the disorder correlation potential in the re- 

1 /3 

stricted regime ((l?oo(0)^)i^ -1) < 1 (seethe discussion 
above Eq. (90)). Next, we must calculate also the replica 
symmetry-breaking solutions of the free energy (45) hav- 
ing more than one discrete step. To solve the minimum 
problem in this case is rather difficult. Therefore, we 
restrict us first to the determination of the continuous 
symmetry-breaking solutions. 



C. Continuous symmetry breaking 

Finally, we look for solutions a{s) of (50) which are 
continuous. In this case, we can solve the stationary 
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equation by a partial integration of B[A(s)] in (49) re- 
sulting in [39] 



ds'a'{s')g'{A[s']). (65) 



Here we assumed that cr(s) = constant for s > Sc- By 
taking two derivatives of (50) we obtain that a{s) fulfills 
the following equation 

a'is) = -2a'(.)a(.)^/^i|§J^-i,,'(A[,s]) . (66) 

Similar as in the case of the one-step symmetry-breaking 
solution we can neglect the longitudinal component in 
f/(A) (47) being a factor cge/cn <C 1 smaller than the 
transverse term in .g(A) (see the discussion below (57)). 
We point out that this is true irrespective of the value of 
A. Using (58) we obtain two solutions of (66) by taking 
once more the derivative with respect to s. This results 
in the following solutions of (66) 



1. cr'(s)=0, 



(67) 



C44 



a{s) s 



A(s) l + iA(s) l + fA(s) 



l + |A(.) + iA2(.) 



(68) 



By inserting (66) into (68) we obtain for the second type 
of solutions 



:(I?(0)yl)i/3 



l + tA(s) 



5/3 



l + |A(s) + iA2(s) 



(69) 



Finally, we determine the constant defined in (65) 
where a{s) = constant for Sc < s < I. By using (50) 
we obtain 



ct(sc) 



V (2B[A(s,)]) 



(70) 



With the help of (68) we obtain 



l + !A(s, 



l + |A(,s,) + iA2(s,) 



xUaaV2(.,) 




which leads with (69) to 



(2?(0)A) 



iAAi/2(s,) l + iA(s,) 



1/2 



l + tA(s,) 



(72) 

Under consideration of (51) we obtain that (72) can be 
solved only for 'D{Q)A > 1. Furthermore, by taking into 



account (69) in the case of A = which is the same equa- 
tion when taking only the most leading A-term for A — > 
in (58), we obtain in this limit no solution of (69). As 
mentioned above, this corresponds to the marginality of 
a string in d = 3 dimensions in an impurity background. 
Due to the non-quadratic polynomial behavior of the ex- 
pressions above, it is not possible to get simple analytic 
solutions in the whole A-range. Therefore, we shall solve 
(67), (69) and (71) for small A <C 1. This corresponds 
to the restriction ((I?(0)A)i/^ - 1)/A < 1. We obtain in 
this range the following solutions: 







for < s < 



Ais)^{ if [(I?(0)A)i/3s_i] for 
.^[(P(0)A)V3_i]2 for 



(D(0)A)i/3 



(15(0)4)1/3 
< S < Sc 



Sr < S < 1. 



with 



(I?(0)A)i/3 



17^2.n^(( ^(0M)^/^-l)^ 
(I?(0)A)7/3 



X>2(0) 



(73) 



(74) 



Finally, we can calculate the free energy A/var for the 
replica symmetry-breaking solution (73) by using (45), 
(58) and (68) for {{V{0)Ay/^ - 1)/A < 1. We obtain 



A/v, 



knT, 



V^iO)\l~{V^{0)A^r'/^] (75) 



identical with the free energy of the one-step replica 
symmetry breaking solution (64). (75) is valid for 
((I?(0)A)i/^ - 1)/A < 1 in the case of a Gaussian dis- 
order potential. But one can generalize the calculation 
above to obtain the validity of (75) in the smaller regime 
((2^oo(0)Aoo)^^'^ — 1) ^ 1 irrespective of the disorder 
correlation function (sec the discussion above Eq. (90)). 
Summarizing, we obtain a saddle point of A/var which is 
symmetric in replica space for 2?oo(0)^oo < 1- In the case 
of I'oo(0)^oo > 1 we obtain a replica-symmetric solution, 
a one-step replica symmetry-breaking solution and also a 
continuous replica symmetry-breaking solution appears. 
To get more insight into the true minimum, we have to 
consider the stability of the various saddle point solutions 
in this case. 



VII. STABILITY OF SOLUTIONS 

In this section we determine whether the various solu- 
tions for the fluid phase discussed in the last section are 
stable in a sense specified below and whether we have 
to take into account also higher-step replica symmetry- 
breaking solutions to get a stable saddle point. A typical 
example of an exactly solvable system with finite-step 
replica symmetry-breaking saddle point solutions which 
are not stable is a string in two dimensions with a S- 
impurity correlation function resulting in an unphysical 
negative variance of the free energy with respect to disor- 
der averaging [51]. This negative variance is vanished in 
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the case of the infinite or continuous repHca symmetry- 
breaking solution. Mezard and Parisi [36] show two dif- 
ferent ways to obtain a theory which includes replica sym- 
metry breaking such as (45)-(50) for random manifolds. 
The first approach consists of a large TV' expansion of 
the partition function where N' is the number of com- 
ponents of the fields. There is only a slight difference 
between the large N' approach and the variational ap- 
proach used above. In the saddle point equation of the 
large A'^'-approach we have to substitute ^{x) in (45)- 
(50) by A{y/x) / (kBT)"^ for the application of this ap- 
proach to the fluid phase of the vortex lattice. This is 
discussed in Appendix B. The large N' expansion con- 
sists effectively in a saddle point approximation in suit- 
able chosen auxiliary fields [36] . The stability of solutions 
of these equations consists in going one step further to 
the quadratic expansion of the action in these auxiliary 
fields with the requirement that the partition function 
calculated from this saddle point approximation is not 
divergent when integrating out the auxiliary fields. It 
was shown by Carlucci et al. [40] that continuous replica 
symmetry-breaking solutions calculated in the last sub- 
section are generally stable in this sense. This is reviewed 
by us in Appendix Bl. Due to the smallness of TV' = 2 in 
the vortex lattice system we do not think that the large 
A'^'-expansion is appropriate in our case. 

We derived (45)-(50) by another way also stated first 
by Mezard and Parisi [36] via the variational approach in 
(36). It is clear that in this case we should require for the 
eigenvalues of the matrix built of the second derivatives 
of AFvar with respect to the self-energies CTq/j that these 
are all positive in the stationary point. Here we take 
further into account the symmetry of and (44) in 
the variation of the free energy. The concrete derivation 
was carried out by Sasik in Ref. 53. Starting from his 
expression for the Hessian we carry out in Appendix B2 
a similar stability analysis as was done in the large TV' 
case by Carlucci et al. in Ref. 40 summarized in Ap- 
pendix Bl. We also find in the variational approach that 
the continuous symmetry-breaking solutions are gener- 
ally stable which means that all eigenvalues of the Hes- 
sian are larger than or equal to zero. Furthermore, we 
show in Appendix B2 that the full Hessian has positive 
or zero eigenvalues if and only if the replicon sector con- 
sists of positive or zero eigenvalues. Thus it is enough to 
consider only the replicon sector for stability. The lowest 
eigenvalues in the replicon sector [40] are given in (B17) 
where k = I ~ r + 1. Here / is replaced by the disorder 
function V in our case and L'j,; is given in (B7). Go is 
the transversal Green with zero self-energy and Ga, 
are the value of the transversal full Green function and 
gap function in the Paris block l<a<i?-|-lina Parisi 
hierarchy of level R. 

We now come to a discussion of the stability of the sad- 
dle point solutions in the fiuid high-temperature phase in 
the symmetric form given in Section VIA and the one- 
step replica symmetric form given in Section VIB. 

First, we consider the saddle point solution given in 



Section VIA for the replica symmetric form. The most 
relevant replicon eigenvalue A(0; 1, 1) is given by (B17) 



3(0) .9'(0) (76) 



A(0; l,l)oc 1 + 4^^4^ i?" 



Here the proportionality factor is positive. The stability 
criterion A(0; 1, 1) > leads to 



^(O)^ < 1 



(77) 



Note that 'D{0)A = 1 corresponds to mi = 1 in the one- 
step replica symmetric solution (61). 

Next, wc consider the stability criterion for solutions 
of the stationarity condition (50) in the one-step replica 
symmetry breaking form. Here, we obtain the lowest 
replicon eigenvalues from (B17) 



A(0;1,1) cx 1 + 4 
,fcsT 



(78) 



X L»" 2 



5(Ai) + — [.g(0)-5(Ai)]M5'(0) 



A(l;2,2) cx 1 



2 5(Ai, 



3'(Ai). 
(79) 



Because g{0) is divergent (58) we obtain for the stability 
criterion A(0; 1, 1) > 



7711 < 1 . 



(80) 



By using (58), (61) and (62) we obtain A(l;2,2) = in 
the leading order in Ai ^ 1,1/^4^. This can be seen 
much easier without using the solutions (61) and (62) 
from equations (59) (60). By taking the square of Eq. 
(59) times the inverse of Eq. (60) we obtain under the 
consideration of D"{x) = 2D'{xY/D{x) which we ob- 
tain from (19) that A(l; 2, 2) = when using (60) in the 
leading order in Ai. This now gives the possibility to 
calculate the stability criterion also in the non-leading 
order in Ai. We obtain 



1 ^ 2 



{ksTfal 1 



A? 



larcsinh 



Al 



1/2^ 



Al 



1 + Ai/4y 



^ [DTi^kBTgiA^ 
D{2kBTg{Ai)/v) 

which results in 

{kBTfal 



1/2 



(81) 



A(l;2,2) cx 1 - 
X D" ( 2^0 



(l-H5Ai/4) 



cUa^ {Ai[l-f Ai/4]}3/2 
(Ai))<0. (82) 



when taking the correlation function (19) into account. 
Summarizing, the one-step replica symmetry-breaking 
solution for the correlation function (19) is unstable. 
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Nevertheless, this instability is very weak for A ^ 1 
which is the interesting region near the glass transition 
line. More generally, we show in Appendix C that all 
finite-step replica symmetry broken solutions of the sad- 
dle point equation (50) are unstable for the Gaussian 
disorder correlation hmction (10). In summary, we have 
shown that the stable self-energy matrix for 'D{Q)A > 1 
has the continuous replica symmetry broken form de- 
rived in Section VIC corresponding to the VG phase. 
For 'D{Q)A < 1 we obtain that the full replica symmet- 
ric solution derived in Section VIA is stable. This phase 
corresponds to the vortex liquid VL. The glass transition 
line between VG and VL is determined by 'D{0)A = 1. 

It is very difficult to solve the saddle point equation 
(50) for a general disorder potential K{q). Neverthe- 
less, the glass transition line I?oo(0)Aoo = 1 where I?oo 
is defined in (18) is valid in the general case. We point 
out that equation (81) leading to the instability of the 
one-step replica symmetry-breaking solution in the case 
of the effective Gaussian disorder potential (10) is still 
valid irrespective of the correlation potential K{q). In 
general, the one-step replica symmetry-breaking solution 
is stable if 

7?i[Ai]ki [a'cxp [~q^kBTg{/^i)/2v\ \ < 1 (83) 



with 



X 4 arcsin 



(1 + 5A/4) 
(A(l + A/4))3/2 

'ai/2\ / A 



(84) 



1/2- 



^1 [/] = 



^2 J \^l + A/4 



(85) 



where we get an unstable one-step replica symmetry- 
breaking solution if and only if 771K1 > 1. Next, we con- 
sider the existence of the continuous replica symmetry- 
breaking solution. For their existence it was crucial that 
in (69) the right hand side was larger than zero. The 
corresponding equation without restrictions on the cor- 
relation function is given by 



s (I?oo(0)^oo)'/' = r;2[A(s)]K2 [Kcxp [-q^B[A{s)]/2] 

(86) 



with 



l + |A(s) 



5/3 



^2 [/] 



22/3 



l + fA(.) + iA2(,s)' 



(87) 



4/3 



To solve this equation for a general correlation function 
K{q) without further approximations is not an easy task. 
Nevertheless, the condition that a continuous replica 
symmetry broken solution exist is given by the possibility 
to solve (86) for Sc resulting in 



miMsc)] 1^2 Kcxp [-q^B[A{s,)]/2] 



> 1. 



(89) 



Quantities like ki,K2 are well known quantities in 
probability theory. The corresponding quantity called 
kurtosis measures in one dimension the curvature of a 
probability function compared to the Gaussian proba- 
bihty function. This is also vahd for (85), (88). The 
Gaussian correlation function (10) has ki, K2 = 1. A cor- 
relation function with a sharper tip and longer tails like 
~ exp[— g"] for a < 2 has ki,K2 > 1- For a > 2 which 
is more flat near the origin with a shorter tail than the 
Gaussian function has ki, K2 < 1. 

Next, we specialize the stability condition (83) and 
(89) to the vicinity of the glass transition line where 
Ai/2^00/2 < 1. By using (61), (62) and (73) one can 
derive the validity of this condition by ((2?oo(0)j4oo)^^'^ — 
1) <C 1 irrespective of the disorder correlation function. 
In this regime, we obtain by an expansion of K around 
the origin which is justified due to A^/^Aao/2 <C 1 



kexp[-q^kBTg{Ai)/2v] 



(fcsr)V(Ai) 



1^2(0) 



KiO) 



f 9 



KiO) 



f 9 



K{0) 



K2 



k{0)exp[-q^B{A{s))] 



i?2(A(.)) 



(90) 



(91) 



i^2(0) 



K{0) 



V9(g2) 



K{0) - 



f 9 



Va(g2) 



K{0) 



One can derive the simple identity 



^i[K]=k{0) 



d 



9(g2) 



K{0) ^ 



d{q'^ 



k{0) (92) 



where is the kurtosis (85) built with the disor- 

der correlation function K (7) in position space. By 
using 7?i[A] « 1 + 20 A/24 and 772 [A] « 1 + 17A/12 
we obtain the simple rules in Table I in the regime 
((X>oo(0)Aoo)^/^ - 1) < 1 for the stable saddle point 
where Aoo has to be taken at the transition point. 
We obtain from Table I a small transition region at 
1 - 20/6A^ < Ki[K] < 1 - 17/6A^ where a higher- 
step replica symmetry broken saddle point solution of 
(36) should give the best free energy. We expect that 
this finite-step replica symmetry-breaking solution leads 
also to a third order glass transition in this range. 

By using A^o — -ic^^iO?^'^^ / a^ksT k. 6/27rc| (see the 
discussion below (63)) we obtain for the high magnetic 
field part 6 ^ 0.5 of the glass transition line the following 
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Kl[K] 


< 1 - 


> 1 - 


saddle point 
order of transition 


one-step breaking 
third order 


continuous breaking 
third order 



TABLE I: Stable saddle points of (36) as a function of the 
kurtosis ki (K) (85) of the disorder correlation function K in 
real space. The second line of the table denotes the character 
of the stable solution of equation (50) . The third line stands 
for the order of the VG-VL transition. 



suTiple result: 

When the kurtosis ki of the positional disorder correla- 
tion function K is smaller than the kurtosis of a Gaussian 
function (flatter tip, shorter tail), the stable saddle point 
solution of (50) is given by a one-step replica symmetry 
broken solution with free energy (64) in the VG phase. 
We obtain a third-order glass transition. When the kur- 
tosis Ki is larger or equal the kurtosis of a Gaussian func- 
tion then we have a continuous replica symmetry broken 
solution with free energy (75) in the VG phase in accor- 
dance with the one-step replica symmetric free energy 
(64). According to Table I we obtain that for lower mag- 
netic fields the border in the disorder function space of 
one-step replica symmetry breaking solutions and conti- 
nous replica symmetry breaking solutions moves to lower 
kurtosises. 



VIII. SOLID PHASE 

In this section we determine the free energy in the solid 
phase. This system corresponding to a string lattice in 
a random potential was discussed in Ref. 39. Here, wc 
reconsider it where we took more emphasis on the de- 
termination of the free energy of the vortex lattice in 
the low-temperature phase than the former work. For 
i^var we use again the approximation (40). For deriv- 
ing this expression one has to consider other arguments 
than in the fluid phase below (40). First, we use that 
the saddle point Green function calculated with (40) ful- 
fills A:bT|Gq^(0) — Gaf}iaei)\/v ^ justified in Ap- 
pendix A. Here ae^ is a nearest neighbor vector in the 
xy-plane. Similar as in the considerations below (14) 
we can restrict us to the diagonal summands x = x' 
corresponding to (40) where non-diagonal terms being a 
factor {Bafj + ^^)^/^/a smaller because as is shown in 
Appendix A Bap ^ for almost the whole range of 
replica indices. We point out that the a, /3-range where 
this inequality is not fulfilled is important for the long 
distance behavior of the lattice fiuctuations beyond the 
random manifold regime [39] (see (103) below). Nev- 
ertheless, due to the a, /3-sum in the various free energy 
terms in (38) one can show that these contributions to the 
free energy are negligible where the inequality BafS <C 
is not fulfilled (see the continuous version (45) of (38) and 
the inequality (A2)). Under the considerations above, wc 
arrive at the disorder Hamiltonian (40) as the basic dis- 
order Hamiltonian for the solid low-temperature phase. 



It is well known and can be also shown by a similar 
analysis as was done for the fluid phase in the last sec- 
tion that finite-step replica symmetry breaking solutions 
are unstable but the continuous replica broken solution 
exist which is stable as is shown by Carlucci et al. in Ref. 
40 and Appendix B. This breaking of the replica symme- 
try corresponds to a glassy phase. We now calculate this 
replica broken solution. The calculation is similar to the 
calculation of the continuous symmetry broken solution 
in the fluid phase carried out in Section VIC. As was 
done before for the fluid phase we can restrict us to the 
transversal fluctuations in the displacement fields u be- 
cause cii 2> C66. Then we have to calculate g{A) where 
we restrict us to the two lowest-order expansion terms in 
A <C 1. By using (12) we obtain 



5(A) 



1 



C44 Vbz jbz 



d kdk^ 



1 

[2-2cos(fc3a3)] + 3^Ej2 



- 2cos(fcia)] 



-A 



^ 0.21- 
C44 \ 16 



(93) 



Here we use the same approximation as in Ref. 32 which 
means o^cqq / C44 w 4/7r. 

First, we determine the two solutions of (50) corre- 
sponding to (67) and (69). This results in 



1. 
2. 



a'(,s) = 0, 
A(s) 



^(I?(0)A)V2.3/2 



(94) 
(95) 



Instead of (71) for Sc m the case of the fluid phase we 
find for the solid phase 



=1/2 



31/2 



p2(o) (I?(0)A)"3/2 . 



(96) 



This value was calculated by using (70) with the approx- 
imation V{2B[A{sc)]) « V{0) valid for c| < 1. This is 
correct in the vicinity of the melting line [32]. Summa- 
rizing, we obtain for A(s) 



A{s) 



^(I?(0)A)l/253/2 foj. s<S^, 

V%0)(V{0)A)-'^ for Se < s < 1 , 



33/ 



(97) 



From yl 3> 1 near the melting line (see the remarks below 
(63)) we obtain that A(s) ^ 1 is in fact fulflUed in the 
magnetic temperature regime, we are interested in (note 
that 'D{0)A < 1 on the melting line for temperatures 
larger or in the vicinity of the glass transition line). Fi- 
nally, we calculate the free energy corresponding to (75) 
by using (45) with (93). For A ^ 1 we obtain 



A/v. 



ViO) 



1 



20 



(98) 



Here we neglect energy terms coming from the first term 
in (45) corresponding to the kinetic part which is a fac- 
tor '--^ 1/10 smaller. When comparing (98) with the free 
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energy of the quadratic disorder case (15) of Section IV 
wc obtain that only the second term in (98) is different. 
This term should cancel the first term in (98) for lower 
temperatures resulting in a vanishing of the reentrant 
behavior of the BG-VG, BG-VL line in the quadratic 
disorder case. 

It can be seen from the derivation above and also Ap- 
pendix A that the actual form of the self-energy func- 
tion A(s) depends on the form of the disorder corre- 
lation function not only by one small parameter. For 
s = Sc we have B{A[sc]) <IC C'^ but for s <C Sc we have 
B{A[sc]) ^ S,'^ which means that the form of the whole 
correlation potential is important when solving the sad- 
dle point equation (50). This makes it difficult to solve 
this equation in general. Nevertheless, for disorder corre- 
lation functions in the vicinity of the effective Gaussian 
disorder correlation function (10) we think that the result 
(98) should not be much changed. 

IX. OBSERVABLE CONSEQUENCES 

Let us now apply the results obtained above to find the 
BG-VG line and the glass transition line of YBCO. The 



entropy and magnetic induction jumps over the transi- 
tion lines will also be discussed. We saw in Section VII 
that the form of the local minimum of the variational 
free energy (36) in the high-temperature phase depends 
not on the kurtosis of the disorder correlation function 
K where the results are summarized in Table I. Here the 
Gaussian disorder correlation function with = 1 

separates the regime where we have a local minimum of 
Fvar of the one-step symmetry-breaking form (Ki[iir] < 1) 
and of the continuous replica symmetry-breaking form 
(Ki[iir] > 1) for large A^o- 



The free energies for both regimes coincide given by 
(64) or (75), respectively. For the solid phase we ob- 
tain the continuous symmetry-breaking solution (98) for 
the Gaussian correlation function. For disorder correla- 
tion functions in the vicinity of the Gaussian correlation 
function we can use this free energy as a first approxima- 
tion for the free energy of a general disorder correlation 
function. Taking into account (45), (46) we obtain 



A/, 



T-*-0 
var 



V{0) 



1 - ^vmvmr 



BG phase , 



l?co(0) [l- (Poo(0)A^)"'/'l e[Poo(0)A^-l] VG-VL phase 



(99) 
(100) 



in the regime near the melting and glass transition line. 
The free energy is given by Fa w iV/var = A^(/var(0) + 
A/var) (45), (46) where the disorder part of the free en- 
ergy A/var is given by (99) in the solid phase and (100) 
in the fluid phase. The intersection criterion correspond- 
ing to (29) in the quadratic approximation in the disorder 
strength which determines the BG-VG, VG-VL line reads 

_ 0^ (1 - b)^ 3.9 ■ 10-5 ^_(2/fc^T)(AfIr°-AfIr-) 

(BG - VG, BG - VL line) . (101) 

Without disorder we have shown in Ref. 32 that the melt- 
ing criterion (31) is equivalent to a Lindemann criterion 
where the Lindemann parameter is given by w 0.18. 
There are many papers which used Lindemann-like cri- 
teria also to determine the disorder induced BG-VG line 
[21-23, 25, 26]. In Ref. 24 Mikitik and Brandt even tried 
to derive a Lindemann-like criterion for the BG-VG, VG- 
VL line from an intersection criterion similar to the one 
used here. Because these Lindemann-like rules do not 
look similar to our microscopically derived melting crite- 
rion (101) we do not try to go further in this direction. 



As derived in section VIB, the glass transition line 
which is the border of the replica symmetric solution of 
the Mezard-Parisi variational calculation and the one- 
step replica symmetry-breaking solution, where the sta- 
bilities was discussed in section VII, is determined by 
nil = 1 in (61) resulting in 

Voo{0)Aoo = 1 (VG - VL line). (102) 

This equation corresponds to the depinning temperature 
Tdp of a one-dimensional string in three dimensions in a 
random environment [1] . The Larkin length Lc is defined 
by the length where we have a coherently pinning of the 
string which means u^{0, Lc) = C'^ ^ £,ab with 

u2(L,L3)= ((u(L,L3)-u(0,0))2). (103) 

When temperature fluctuations become larger there is a 
softening of the impurity potential which is important 
for the length of the coherently pinned vortices. This 
correction is important when this fluctuation length be- 
comes equal to Lc calculated for T = 0. This depinning 
temperature is given by (102) [1]. We mention that it is 
difhcult to distinguish experimentally by diffraction ex- 
periments in which of the two classes ki [/f] > 1 — 17/6A^ 
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FIG. 4: The BG-VG, BG-VL first-order transition lines 
B,n{T) given in (101) for JTc-pinning (upper figure) and Sl- 
piiniirig (lower figure). The solid lines are calculated with 
parameters for do and ^' which gives one of the best fits to 
the experimentally determined [6] BG-VL line (square points) 
within the pinning mechanism (27rdo Cab/C' — 1-5 • 10~^ and 
Cab/C' = 1-59 for (JTc-pinning, 27rdoCab/^'^ = 1-32 ■ 10"'' and 
^ab/C' = 1-49 for (5/-pinning) . Dotted curves are calculated 
by a variation of these parameters given by = (1 ± l/2)do 
and 1/5'* = (1 ± 1/2)^/^/5' where only one parameter was 
varied. We wrote that parameter at the curve. The vertical 
markers denote the intersection points of the glass transition 
line and the BG-VG, BG-VL line named GP. From Clausius- 
Clapeyron equation (105) the critical points CPs are deter- 
mined by an extremum of the BG-VG, BG-VL lines 



or < 1 — 20/6A^ the disorder correlation poten- 

tial _?C of a given experiment belongs. In both regimes 
we obtain u^(0, L3) cx {kBT)L^/ diO?' in the VG-VL 
phase (the proportionahty constant is different for both 
regimes, see also the discussion in Appendix A). This is 
reasoned in the vanishing support of (j{s) in the vicinity 
of the origin, fact for the the one-step rephca symmetry- 
breaking regime ni[K] < 1 — 20/6 (56) as well as the 
continuous rephca symmetry-breaking regime ni [K] > 
1 - 17/6A^ (68) and (73). It means that thermal fluc- 
tuations are dominant over disorder fluctuations. Note 
that u^{L,L3) for L ^ diverges in the VG as well as 
the VL phase characteristic for defect dominated phases 
also seen before for the system without disorder. For 
the BG-phase in the random manifold regime we obtain 
u^{L, L3) (X (L^-|-L|)^/^ in accordance with former calcu- 
lations [39] . The derivation beyond the random manifold 
regime where the lattice structure is important leads to 



1 ' \ ' r 

5Z-pinning 




TIK] 



FIG. 5: Phase diagram for YBCO. Solid lines represent the 
theoretical determined phase transition lines between the var- 
ious phases calculated for (5Z-pinning with 27rdo 5ab/C' ~ 
1.32 • 10"® and ^ab/£.' = 1.49 corresponding to the solid line 
in the lower picture in Fig. 4. The glass transition line VG- 
VL was calculated from (102). Square points represent the 
experimentally determined phase diagram of Bouquet et al. 
[6]. 



?/2(L,L3)cxlog(L2+L2) [39]. 

Mikitik and Brandt found in Ref. 23, 24 that their ana- 
lytical derived BG-VG, BG-VL curve is a function of the 
Ginzburg number Gi = 32TT\Xab{0)^c{0)T,/cj)lU{0)? , 
b, T/Tc, and the disorder function 'D{0). This can be 
also shown easily for the melting curve (101). Note that 
the disorder constant D in Ref. 23, 24 is a function of b, 
T/Tc, Gi and our disorder function I?(0). 

In Fig. 4 we show the BG-VG, BG-VL curves given by 
(101) for various values do, The upper curves are cal- 
culated with a (5Tc-pinning correlation function (8), the 
lower curves for a (5/-pinning impurity correlation func- 
tion (9). For clearance we do not show the critical points 
CP on the melting line in the figure which are charac- 
terized by zero entropy jumps AS'; per double layer and 
vortex over the transition line. These points can be easily 
marked in the figure since they correspond to the extrema 
of the melting line i?,„ due to the Clausius-Clapeyron 
equation given in (105) below. 

The intersection point of the glass transition line BG- 
VL which is calculated by (102) with the BG-VG, BG-VL 
line is denoted by GP in the figure. The square points 
with the dashed line denotes the experimentally deter- 
mined BG-VG, BG-VL line of Bouquet et al. [6] shown 
also in Fig. 2 for comparison. In the 5Tc part of the figure, 
we find no solutions for the equation (101) near T « T^. 
The parameters of the straight lines in the figure are cho- 
sen in such a way that we reproduce in one of the best 
ways the form of the experimentally melting line of Bou- 
quet et al. [6] and also the position of the experimentally 
found CP and GP. These experimentally chosen param- 
eters arc 2Trdo^lJC'^ = 1.32 • lO"'' and ^ab/C = 1-49 for 
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5Z-pinning, 2^do ^^^/C'' = 1.5 • 10"^ and UK' = 1-59 
for 5Tc-pinning. Thus, wc obtain that the correlation 
length S^' of the disorder potential almost corresponds to 
the coherence length S^ab of the superconductor. The rea- 
son that U /C is larger than one could be due to lattice 
influences on the effective broadening of the vortex (see 
the notes below (9)). Finally, we mention the similar- 
ity of the do parameter values in the Parisi case and the 
corresponding values in the quadratic disorder case of 
Section IV. 

The curves of representative variations of these almost 
optimal parameter values arc shown by the dotted curves. 
We obtain from the figure as was also the case in the 
second-order perturbativc discussion in Section IV that 
the JTc-pinning curves fits less to the experiment than the 
(5Z-pinning curves. This comes mainly from the smooth- 
ness of the disorder parameter d{T) in (9) as a function 
of T resulting in the slow variation of the transition line 
Bm{T) seen in the upper part of Fig. 4. From Fig. 4 we 
obtain that the glass intersection point GP and the crit- 
ical point CP does in generally not coincide. This was 
just mentioned in Ref. 13 for BSCCO where in the exper- 
iments this difference is not seen yet maybe because of 
experimental uncertainties. 

One of the most interesting results of our calculation 
is that the reentrant behavior of the melting line and 
the experimentally not seen low B parts of the BG-VG 
curves in the quadratic disorder calculation of section 
IV (see Figs. 2 and 3) vanished in the Parisi approach. 
It is remarkably that the large descend of the curves 
in the direction to lower temperatures in Fig. 2 within 
the quadratic approach arc smoothed within the Mczard- 
Parisi approach such that the BG-VG transition curves is 
almost horizontal. There arc various forms of the BG-VG 
lines in the literature. One of the reasons of the differ- 
ences in the various experiments comes presumably from 
the strong dependence of the BG-VG line on the dcpin- 
ning function d{T) via an exponential behavior in (101). 
Any pcrturbational effects like surface effects or twin- 
ning areas in the crystal can change the functional form 
of the curve at small temperatures easily. We note that 
especially the strong dependence of the BG-VG curve on 
small variations of the disorder correlation length ^' hav- 
ing its reason in the quadratic dependence of ^' in A (63) 
which is contained as a third-order summand in the free 
energy of the solid phase in (98). The form of the free 
energy in the solid phase is the most dominant factor for 
the form of the BG-VG curve in the vicinity of the criti- 
cal point CP. This is in contrast to the free energy in the 
high-temperature phase which just grows in importance 
beyond the glass intersection point GP but is still small 
compared to the free energy part of the solid phase. 

In Fig. 5 we show the whole phase diagram for the pa- 
rameters of the solid curve in the lower ^/-pinning part 
of Fig. 4 (solid curve). For comparison (square points 
with dashed curve) we show also the experimentally de- 
termined phase diagram of Bouquet et al. of Ref. 6. Both 
phase diagrams look rather similar except that the CP 



and GP points of the theoretical determined phase dia- 
gram lies a little bit lower in temperature in comparison 
to the experimental ones. The upper curve between the 
VG-VL phases show the glass transition line calculated 
by (102). There is a small discrepancy in the slope of the 
line between theory and experiment. Note that for A w A 
which is the case for BSCCO [41] we get that V^{Q)A^ 
is in fact independent of the magnetic field B resulting in 
a vertical glass transition line. This is in good accordance 
to the experiments [13]. 

Next, we calculate the entropy and magnetic induction 
jumps over the BG-VG, BG-VL first-order line. Denoting 
the spacing between the Cu02 double layers by we 
obtain for the entropy jump per double layer and vortex 
over the BG-VG, BG-VL line 

« UbTJ-^j^ HZl-^/Zl^^] (104) 

and a corresponding equation for the glass transition line. 
Now we make use of the Clausius-Clapeyron equation 
which relates the jump of the entropy density of a first- 
order transition line to the jump of the magnetic induc- 
tion by 

Here iJ,„ is the external magnetic field on the BG-VG, 
BG-VL hue. Because B - iJc2(T) for YBCO we can use 
H ^ B in the Clausius-Clapeyron equation (105). Equa- 
tion (105) is not appropriate for a numerical evaluation 
because of the vanishing of the denominator at the sad- 
dle points of the BG-VG, BG-VL line which are canceled 
due to zero points in the numerator. By using the inter- 
section criterion for the transition line we can transform 
(105) to 

AB « fesT™ — -|- ln[Z^^°°/Z^^"] . (106) 

V OBm 

This equation can be also derived from thermodynamical 
relations under the considerations AB / Bm "C 1 which we 
also used by taking the intersection criterion for the free 
energy and not for the corresponding Gibb's potential in 
this paper [54]. 

In Fig. 6 we show AS"; and AB for the parameters 
used in Fig. 5 over both lines. We show in the upper part 
of the figure ASi with experimental points of various 
torque and SQUID experiments (circles [55], squares 
[56], triangles [6]) for the entropy jump over the BG-VG, 
BG-VL line. For the parameters used in Fig. 5 we obtain 
a value for the CP of &%K . In the lower part of Fig. 6 
we show the magnetic induction jumps AB. The square 
and circle points are experiments (circles [55], squares 
[57]). Finally, we note that our curves of the entropy 
and magnetic induction jumps over the BG-VG, BG-VL 
line is in qualitative agreement with similar curves 
calculated within the Ginzburg-Landau approach for 
YBCO by Li and Rosenstein [20]. The main difference is 
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FIG. 6: In the upper figure we show the entropy jump AS; 
per double layer and vortex according to (104). The points 
in the figure are entropy jumps determined by experiments 
(circles [55], squares [56], triangles [6]). In the lower figure 
we show the magnetic induction jumps A_B calculated by the 
help of (106). Experimental points in this figure are from 
Ref. 55 (circles) and Ref. 57 (squares). The solid curves 
in both figures correspond to the jumps over the BG-VG, 
BG-VL line. We used for the whole figure parameter values 
27rdo Cfb/C'^ = 1.32-10"'' and (,abK = 1-49 in correspondence 
to the parameter values in Fig. 5. 



that they obtain a zero point in the magnetic induction 
jump curve in the vicinity of the critical point which has 
its reason in the reentrant behavior of their calculated 
melting line to second order in the disorder potential 
(see (105) by taking into account that dBm/dT is infinite 
at the reentrant points). We expect, as was also the case 
in the elasticity approach used here, that this zero point 
vanishes when going beyond second-order perturbation 
theory leading to the vanishing of the reentrant behavior. 



Finally, we come back to a discussion of the scenarios 
of the phase diagram for YBCO given in the introduc- 
tion of this paper. We did not find a slush phase within 
our numerical examinations of equation (101) during this 
work irrespective of the parameter range. This in ac- 
cordance to the Ginzburg-Landau calculations of Li and 
Rosenstein in Ref. 20. This means that our phase di- 
agram is only in accordance with the second scenario 
of a unified BG-VG, BG-VL first-order line discussed 
in the introduction of this paper. Due to the contro- 
versy of this phase we cannot determine within our the- 
oretical approach whether it is in fact existent or not. 
It was claimed in the experimental paper [12] that the 



slush phase only exists within a really small doping re- 
gion where the entropy jumps over the first-order line 
between the slush phase VS and the vortex liquid VL 
is two orders smaller than the entropy jumps over the 
BG-VL melting line. The intersection criterion of the 
high and low-tcmpcrature free energy used in this paper 
by a perturbativc calculation in both phases uses the as- 
sumption that the slope difference corresponding to the 
entropy jumps is not too small. This could be the reason 
that we do not see the slush phase. We point out that 
to our knowledge there exist no theoretical model which 
shows without doubt the existence of this phase. 



One of the main findings in this work for vortex model 
(2) is that the order of the glass transition VG-VL is 
of third order irrespective of the form of the disorder 
correlation potential. We point out that a third order 
phase transition having a smooth heat capacity should 
show scaling behavior with a non-trivial fix point in a 
rcnormalization group calculation. Prominent examples 
of third order phase transitions is the non-interacting ho- 
mogeneous three dimensional Bose gas across the Bose 
Einstein transition [58] or the large A^'-limit of the two 
dimensional U{N') lattice gauge theory with a variation 
in the coupling constant [59]. It was noticed in Ref. [58] 
that the heat capacity curves for BSCCO over the super- 
conducting transition without magnetic field looks rather 
similar to the heat capacity curves of the homogeneous 
Bose gas. For YBCO this transition looks more like the 
A-transition of ^He. A discussion of scaling relations in 
higher order phase transitions and their classification due 
to Ehrenfest can be found in Ref. 60. 



Fisher et al. proposed in Ref. 9 a scaling behavior 
of the VG-VL glass transition where they introduce a 
disorder phase correlation length with scaling S^q ^ 
|r — TgI"" in the fluid phase near glass transition tem- 
perature Tg on the VG-VL transition line. This scaling 
proposal was later on approved experimentally via mea- 
surements of the current voltage characteristics over the 
transition region [10]. There are now a number of ex- 
periments [10, 61] and computer simulations [19, 49, 62] 
of various models for superconductors showing also this 
scaling behavior where in most cases the disorder phase 
correlation exponent lies in between 0.8 < ly < 1.7. One 
can connect the phase correlation scaling exponent v with 
the heat capacity exponent a defined by C ~ |T — Tg|~" 
where C is the heat capacity via the hyperscaling relation 
1^(1 = 2 — a. Here d is the dimension of the system which 
means d = 3 in our case. Thus, most of the experimen- 
tally determined and computer simulated systems have 
an a-exponent lying in between —3.1 < a < —0.4. This 
corresponds to a phase transition of order three or even 
higher within the Ehrenfest definition of phase transi- 
tions [60]. 
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X. SUMMARY 

In this paper, we have derived the phase diagram for 
superconductors having their phase transition hnes at 
high magnetic fields near Hc2, such as YBCO. The aim 
was to obtain a unified analytic theory for the BG-VG, 
BG-VL transition as well as for the glass transition lines. 
The model consists of the elastic degrees of freedom of 
the vortices with additional defect fields describing in the 
most simple way the defect degrees of freedom of the vor- 
tex lattice. For the impurity potential we restricted us 
to weak pinning 5Tc and (5Z-correlated impurities [1]. 

First, we have derived the effective low- (12) and high- 
temperature Hamiltonians (13) without disorder in Sec- 
tion III. The low-temperature Hamiltonian consist of the 
well-known elastic Hamiltonian of a vortex lattice where 
defects are frozen out. At high-temperatures, the stress 
fields are frozen out leading to the high-temperature 
Hamiltonian (13). In Section IV we have carried out 
the disorder averaging to second-order perturbation the- 
ory with these low- and high-temperature Hamiltonians 
to find the BG-VG, BG-VL transition line by the ap- 
plication of the intersection criterion. The result given 
in Eq. (29) and displayed in Fig. 2 shows a reentrant 
behavior. The low-_B behavior of the calculated transi- 
tion line was not in agreement with experiment. This 
led us to calculate the free energy in the low- and high- 
temperature phases using the non-perturbative approach 
of Mezard-Parisi. In Section VI we calculated the varia- 
tional free energy in the high-temperature liquid phase. 
We obtain a glass transition from a replica symmetric 
solution corresponding to the vortex liquid VL to a sym- 
metry broken solution corresponding to the vortex glass 
phase VG. The position of the glass transition line ful- 
fills equation (102) describing the depinning transition 
of a string with stiffness 0440^ in three dimensions. The 
degree of replica symmetry breaking of the variational 
Hamiltonian depends on the form of the disorder corre- 
lation function. The high-temperature part of the free 
energy is given by (100). For high magnetic fields near 
Hc2 we got the following result: We obtain a one-step 
replica symmetry-breaking solution when the kurtosis ki 
of the disorder correlation function in position space de- 
fined in (85) is smaller than one. In the case that the 
kurtosis is larger than or equal to one we obtain a full 
replica symmetry broken solution. In both cases we ob- 
tain a third order glass transition line. The Gaussian 
correlation function is the border in the disorder corre- 
lation function space with ki = 1. Corrections to this 
simple rule relevant for lower magnetic fields are given in 
Table I. 

In Section VIII, we calculate the free energy of the 
vortex system in the low-temperature solid phase (BG), 
given by (99). The stationary solution for the self-energy 
matrix in replica space is continuous replica symmetry 
broken. By using the intersection criterion for the low- 
and high-temperature free energies, we calculate the ex- 
pression for the unified BG-VG, BG-VL line given by 



(101). In Fig. 4 we show the unified BG-VG, BG-VL 
line for various parameters for both pinning mechanisms. 
We obtain that ^/-pinning fits much better to the exper- 
iments than (5Tc-pinning. It is seen that the reentrant 
behavior of the second-order perturbation theory carried 
out in Section HI vanished in this non-perturbative ap- 
proach. In Fig. 5, we show the theoretical determined 
phase diagram for YBCO. Fig. 6 shows the entropy jumps 
and magnetic field jumps over the BG-VG, BG-VL tran- 
sition line. Finally, we calculated heat capacity scaling 
exponents a from disorder phase correlation exponents v 
determined from experiments and computer simulations 
via the hypercaling relation across the glass transition 
line VG-VL which is only consistent with a third or even 
higher order VG-VL phase transition line. 

APPENDIX A: JUSTIFICATIONS FOR 
APPROXIMATION OF DISORDER 
HAMILTONIAN (40) 

We restrict us here to the case of transversal fluctua- 
tion where the generalization to arbritrary fluctuations 
is straight forward. That (40) is valid for the high- 
temperature fluid phase was shown below (43). 

In the solid phase, we first have to show that 
A:bT2 |Gf (aci) — GJ(0)|/w <C in the interesting regime 
near the melting line where is the full Green function. 
aei is a nearest neighbor vector in the xy-plane and s is 
a continuous Parisi index. It follows from [36] 

Gf (x) = / Skdk^ e'^-^ Gl (Al) 
Vbz Jbz 

'1 A(s) r ds A(s) 

I(Gj)-i + A(s)+X ^(Gj)-i + A(.). 

that the nearest neighbor fluctuations are in fact much 
smaller than the nearest neighbor distance a when taking 
into account (93), (97), c| < 1 and D{A)A < 1 near 
the melting and glass line which is the regime we are 
interested in. 

Finally, we have to show that B{A{s)) ^ for al- 
most aU s > 0. From (49), (97) with (93), c| < 1 and 
D{0)A < 1 we obtain B{A{s)) < for 

s > 4 . (A2) 

which is almost the whole s-region. The extreme small 
range s ^ c| has no relevance for the free energy result. 
As mentioned above, this small s-region of ct(s) becomes 
relevant only when calculating disorder fluctuations (103) 
beyond the random manifold regime which corresponds 
to distances L, L3 where v?'{L,L^) S> [39]. 

APPENDIX B: STABILITY OF MEZARD-PARISI 
SOLUTIONS 

In this section, we consider the stability criterion of 
the Mezard-Parisi theory in the large iV'-limit and in the 
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Bogoliubov variational method. First, we reconsider the 
derivations of Carlucci et al. [40] for the stabihty condi- 
tions in the case of the large A'''-liniit. Then wc derive 
the corresponding stability criteria in the variational ap- 
proach considered in Section VI. To our knowledge this 
was not done before in the literature. 

In order to compare the vortex lattice theory with two 
component displacement fields with the N' = 2 isotropic 
random manifold theory of Mezard and Parisi we re- 
strict us in the following first to the transversal displace- 
ment fields justified above as a good approximation in 
both phases. A generalization to the full fluctuations is 
straight forward. The difference of the stationary and 
stability expressions in both phases for the vortex lat- 
tice and the isotropic N' = 2 random manifold theory of 
Mezard and Parisi comes then mainly due to a difference 
in the kinetic part of the Hamiltonian Gq ^ in Eq. (39). 
As was shown in Ref. 36 the saddle point equation of both 
approaches looks rather similar except that in the large 
A'''-limit the saddle point equation (50), I?o(x) should 
be replaced by f{x) where f{x) = A{^/x) / (K bT)-^ (we 
take the reversed sign to the Mezard-Parisi definition). 
A(y^) is the impurity correlation function in (32). To 
derive this, Mezard and Parisi insert in the action of the 
isotropic random manifold system, auxiliary fields. By 
integrating out the fluctuating displacement fields of the 
random manifold the large A^'-limit corresponds to a sad- 
dle point approximation in the auxiliary fields. This re- 
sults in the saddle point equation (32). By the definition 
of 

fix) = / da e-"f{ax) (Bl) 
Jo 

Mezard and Parisi obtain (50) for the general variational 
approach where the disorder function Voix) is replaced 
by Rx). 

1. stability in the large TV'-limit approach of 
Mezard and Parisi 

The stability of the stationary solution (50) comes from 
the stability of the saddle point approximation of the 
action in the auxiliary fields. This results in the stability 
matrix [36] (we take into account only the less stable part 
of the stability matrix corresponding to zero moments) 

with 

-viJ ~ ^ ^ ^ / d^kdkz (Ga'j — GaS — G + Gps) 

V Vbz Jbz 

(B3) 

^al -yS — ^ -y^r ^ f d^kdkj (G^^ — GaS — G p-^ + Gps)'^ 

w Vbz Jbz 

(B4) 



and a < P and 7 < i5. Here Gap^k) stands for the 
transversal component of the Green function (39) in the 
case of the vortex lattice or the corresponding Green 
function in the case of the isotropic two component 
N' = 2 random manifold system [36]. The stability of 
the saddle point of (38) fulfilling the discrete version of 
the self-energy equation (50) is given when all eigenvalues 
of the stability matrix (B2) are positive. 

j^aP-nS jg fQ^j. index ultrametric matrix [37]. It was 
shown by Kondor et al. [63] and later on by Temesvari 
et al. [52] that one can divide the eigenvalues of ma- 
trices of the form (B2) in three classes. The first two 
families consist of vectors in the longitudinal sector of 
dimension R + 1 and R anomalous sectors of dimension 
R + 1 depending explicitly on the form of the ultramet- 
ric matrix. Here, we denote R by the level of hierarchy 
of the self-energy matrix cTap fulfilling the stationarity 
condition (50). This means i? = for the replica sym- 
metric solution calculated in section VIA and i? = 1 for 
the one-step solution given in section VIB for the fluid 
phase. 

There is no closed form in the literature for the eigen- 
vectors and eigenvalues of the matrix in (B2) for the first 
two families. Nevertheless, it is able to block diagonalize 
the matrix M given by (B2) in the various sectors [52]. 
Following Temesvari et al. [52] , we denote the size of the 
Parisi blocks as pr, r = 1 . . . R^ where R is the maximum 
level of replica symmetry breaking. We denote po = n 
and = Ij the latter being the size of diagonal el- 

ements. The matrix elements ctq/j, that belong to the 
rth level of replica symmetry breaking are all equal to a 
number denoted by (7^, r = 0, . . . , i?. The replica overlap 
function is defined hy a C\ j3 = r when Uap = Or ■ 

Denoting with < r < i? the basis vectors in the 
first two families. For k = which is the longitudinal 
sector we obtain for the i? -t- 1 basis vectors [52] 

, On / 1 for a n /? r , , . 

("'■^"^^lo for an/3^r. 

The basis vectors for fc ^ corresponding to the 
anomalous sector can be found in Ref. 52. 

The third family of eigenvectors of ultrametric matri- 
ces as for example (B2) is named the replicon sector. It 
consist on several one-dimensional subfamilies labeled by 
r = 0, . . . , R and k,l = r + 1, . . . , R+1. The correspond- 
ing one-dimensional subspaces are eigenspaces with the 
eigenvalues denoted by A(r; k, I). The eigenvectors corre- 
sponding to the basis vectors in this sector can be found 
in Ref. 52. Note that these eigenvectors do not depend 
on the entries of the ultrametric matrix. The eigenvalues 
A(r; k, I) can be generally expressed via the matrix ele- 
ments of the ultrametric matrix [52] . In the case of the 
concrete ultrametric matrix M (B2) one finds [40] 

Hr;k,l) = r,f,,fr)fr^ ^/^Z \\ ^ ^kl (^6) 

2J"{2(KBT/v){gR+i ~ gr)) 
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with 



(Pkdk^ 

BZ 



(B7) 



[(Go 



Here gk corresponds to the value of the transversal com- 
ponent of the Green function Gap (39) integrated over 
the momenta as in (47) or the corresponding random 
manifold Green function with a H [3 = k. The eigenval- 
ues A(r; k,l) ior k = I = r + 1 are the most singular ones 
for definite r. One can show easily [40] that these most 
singular eigenvalues are zero in the case of continuous 
symmetry-breaking solutions as in section VIC for the 
fluid phase. We note that for the stability matrix sector 
of moments unequal to zero there are only eigenvalues 
larger than zero [40] . 

Finally, we sketch the proof given in Ref . 40 that eigen- 
values of the first two families, which is the longitudi- 
nal sector and the anomalous one, has only eigenvalues 
which are larger or equal to the replicon eigenvalues given 
above. With the definitions 



\{Pr -Pr+l) 
Pr - Pr+l 



for r < fc — 1 
for r ~ k — 1 
for r > k — 1 . 



and 



Afe(r) = 



we obtain 



A(r; k, r + 
A(r; r + 1, r 



1) 



1) 



for 
for 



k > r 
r > k - 



det 



Al) =n[Ak(r)-A]det 



r=0 



I + M' 



(k) 



with 



A/' 



AT 



2(Afc(s) ~ A) 



(B8) 



(B9) 



(BIO) 



(Bll) 



where we denote M'^'^'' by the matrix AI'l-l^ and KJl^ is 
a generalized discrete Fourier transform of the ultramet- 
ric matrix M^'') [40]. The matrix Kl" is denoted as the 
kernel for the ultrametric matrix Af"l^'^^ which means 
that M"'^-'^'' is given by M^^s^ = Afc(r) + A:[^A^/2 in the 
longitudinal or anomalous sector k and r, s runs over the 
basis vectors in the k sector. 

We point out that KJ^," can be written as 



/qO' = 4Bfc(max(r, s)) 



(B12) 



where the function Bk can be expressed explicitly by the 
Green functions Gg [40] so does not depend on the dis- 
order function /. One finds 

Bk{r) < and Bk{r + 1) - Bk{r) > . (B13) 



Denoting det^^^/ [M'^"^^] by the determinant of the sub- 
ik) 

matrix of M with lines in S G {0, . . . , i?} and 
columns in S' G {0, . . . , R} where we suppose that 
S and 5' has the same number of elements denoted by 
#5' = #5'. We obtain 



det (m^'^) - AI 



n[Ak(r)-A]^det5,5[ 



T'(k)i 



1-0 



(B14) 

In [40] it is shown that the right hand side is larger than 
zero for A < Minr[Afe(r)]. This means that the eigen- 
values in the non-replicon sectors are always larger than 
the smallest eigenvalue in the replicon sector. We now 
give a more general proof of this fact useful in the next 
subsection: 

This is true if det^,^. [M'^''^] > for S,S' G 
{0, . . . ,R}. We suppose the ordering Sj < Si+i for s e 5 
and similar for 5'. By subtracting appropriate line and 
columns of the matrix where Kj. denotes the matrix 
K^"^ we obtain a matrix where its determinant is built 
purely from its diagonal elements given by 

det5,5'[Kfe] = 4*^Bkimax{s^s,s'#s,)) 

#5-1 

X Yl [5fc(max(sj,s-)) - Sfc(max(.Si+i,s'j_|_i))] . 

8=1 

(B15) 

By using (Bll) with A^' < (B8) for n ^ and (B13) 

we obtain dets,^' [M'*'"'] > for A < Min,.[Afc(r)]. 

It is clear from the considerations above that the var- 
ious sectors especially the longitudinal sector depends 
on the concrete hierarchical structure we choose. This 
means, that we can also get other eigenvalues for the 
various sectors by starting from a given minimal level of 
hierarchy by an appropriate artificial division of the vari- 
ous sectors leading to a larger level of hierarchy. It is clear 
that nevertheless the lowest eigenvalue of the stability 
matrix being in the replicon sector didn't change. Now 
suppose, we try to restrict the stability matrix M'^^''^^ 
to the fc = longitudinal sector of a suitable subdivided 
hierarchy, corresponding to a search of the minimum of 
-Fvar (36) in the self-energy matrices (JafS which are con- 
tained in the Parisi-algebra. 

For a subdivision of blocks we obtain that Bk (r) given 
explicitly in [40] is constant on two blocks in the subdi- 
vided hierarchy originating from the same blocks fc and 
r in the precursor hierarchy. Furthermore, we get dou- 
blings in the eigenvalues A*"' (r) corresponding to the sub- 
division. But this results to det5.5[M^''''] ^ only if 
S does not contain two blocks in the subdivided hierar- 
chy originating from the same block. Then we immedi- 
ately obtain from (Bll), (B12), (B14) that we can always 
subdivide the hierarchy in such a way that the lowest 
eigenvalue in the fc = longitudinal sector is given by 
the minimum of the eigenvalues in the replicon sector 
Minr[Ao(r)]. This means that by restricting the stability 
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matrix M"'^''^^ to the subspace of symmetric self-energy 
matrices in the Parisi-algebra with the constraint (44) we 
obtain that the lowest eigenvalue of the restricted matrix 
]\/[ai3,iS gqyal to the lowest eigenvalue in the replicon 
sector. 



Stability in the variational approach of Mezard 
and Parisi 



In this subsection, we carry out a similar analysis for 
the Bogoliubov variational approach of the Mezard-Parisi 
theory, outlined in Section IV, as was done for the large 
A'''-limit theory in the last subsection. The self-energy 
within this approach is calculated by searching for the 
stationary points of the variational free energy (35). We 
get a stability matrix of this stationary point by taking 
the second derivative of i^var with respect to the self- 
energy matrix under the constraint (44). This was cal- 
culated in Ref. 53. We obtain 



M 



KbT d^F^, 



1 



-L 



(2) 



(2) 



(B16) 



(2) 



This matrix corresponds to the matrix Af^^'^^ (B2) in 
the large A'''-limit approach. Because M'^^'^^ is a ultra- 
metric matrix we obtain by using the rather general con- 
sideration for eigenvalues in the replicon sector of these 
types of matrices [52] 



,i^--n2(KBT/v)[gn+i 



kl 



(B17) 

By comparing (B17) with (B6) we obtain also in the 
variational approach that the most divergent eigenval- 
ues X{r;r -|- l,r -f- 1) are zero in the continuous replica 
symmetry-breaking solutions as was also the case in the 
large A'^'-limit approach. 

In the following we show that in the variational ap- 
proach the eigenvalues of the longitudinal and anomalous 
sectors are larger than zero. We first define the reduced 
stability matrix 



Mred - (L(^))"^M(L(^')-^ 



(B18) 

Here L^^^ is the matrix L^'^p ^g. By using that the ker- 
nel of L^^-* is given by —K^'^ we obtain for the kernel of 
(L('))-i [40] 

= iAl)-'/'m%\Alr'/' KkiAlr' (bi9) 

with 



^(a2)-^/%a^(a°)-V2 



where {A^)rs = A"(r) 5rs and A".(7') is given by 



A°.(r) ^ 



k,r+l 



for k > r + 1 



^r+l.r + l fO'' r > k — 1 . 



(B20) 



(B21) 



We mention that L^^-* has positive eigenvalues which can 
be seen from the positivity of G~^. This has to be as- 
sumed for the stability of -Ftriai, leading to the positivity 
of L^^-* for n because only the first term in (B16) is 
unequal to zero for / = where the second term in (38) 
does not contribute to the stability matrix [53]. Now we 
use that L^^-* is given by A° — Kj,A'^/2 which means that 
has only positive eigenvalues. By using the same con- 
siderations for the eigenvalue equation det[N5J — A] = 
as was done at the end of the last subsection we obtain 
further that all eigenvalues of N° are lower than one. 
This leads to the fact that the denominator (N^)^^ in 
F° (B19) can be expanded in a geometric series. 

The eigenvalue equation for Mj.gj is given by the right 
hand side of Eq. (BIO) with Kl' in (Bll) is substi- 
tuted by the expanded form of {F°Y' (B19). Afc(r) is 
built of the replicon eigenvalues of Mrcd corresponding 
to (B9). By carrying out the calculation of the result- 
ing sub-determinants of sums and products of matrices 
by standard rules (Cauchy-Binet formula) we obtain as 
in the last subsection that the non-replicon eigenvalues 
of Mrcd larger than the smallest replicon eigenvalue 
given in (B17). 

Furthermore, we obtain also with a similar proof as 
in subsection 1 that the projected matrix M"J^'^^ to the 
space of the symmetric self-energy matrices (Tq/3 of the 
Parisi form with the constraint (44) contains the smallest 
replicon eigenvalue. 

Up to now, we have only shown that the results of the 
large A^' approach considered in the last subsection are 
also valid for the reduced stability matrix M^J^''^^ ■ It 
is not clear whether this is also valid for the full stabil- 
ity matrix M^^'f^ (B16) of the Mezard-Parisi variational 
approach. Nevertheless, one normally does not need the 
results above in their general form for a stability analysis 
of saddle point solutions of (36). It is enough for this 
analysis to know the results concerning the positivity of 
the eigenvalues. This can be immediately reached by us- 
ing the defining equation (BIB) of the reduced stability 
matrix and the general conclusions above. 

This leads to the following results for the stabilities in 
the large A^' and the variational approach of the Mezard- 
Parisi theory: 

1. The eigenvalues in the replicon sector are given by 
A(r; fc, I) in (B6) for the large A^'-approach and by 
A(r; k, I) in (B17) for the variational approach. In 
the case that the eigenvalues in the replicon sector 
are all positive in the large A^'-approach or the vari- 
ational approach wc obtain also that all eigenvalues 
of the full stabihty matrices M°'f^''^^ or M"/^'')"', re- 
spectively, are larger than zero. This leads to the 
stability of the corresponding saddle point solution. 

2. The eigenvalues of the continuous symmetry- 
breaking solution are larger than or equal to zero. 

3. The eigenvalues of the stability matrix projected 
on the subspace of variations in the symmetric self- 
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energy matrices aap in the Parisi algebra with the 
constraint (44) are larger than or equal to zero if 
and only if the eigenvalues of the full stability ma- 
trix not restricted to variations in the Parisi algebra 
in both approaches are larger than or equal to zero. 



In (CI) we used that Aq = and mjj+i = 1. The R 
stationarity conditions i9A/var/9Ai = for i = 1 . . . i? 
lead to 



We further note that the eigenvalues A(r; k, I) (B6) and 
A(r; k, I) (B17) of both approaches are proportional to 
each other with a positive proportional constant when 
neglecting the distinction in the effective disorder func- 
tions / and / related by (Bl). 



APPENDIX C: INSTABILITY OF FINITE-STEP 
REPLICA SYMMETRY-BREAKING SOLUTIONS 
IN THE FLUID PHASE 



i=l 



m.1 



.2M^i?'(2i?[A™J). (C4) 



for I = 1 . . . R corresponding to (59) in the case of the 
one-step replica symmetry-breaking solution. The saddle 
point conditions SA/var/^m; = for i = 1 . . . i? lead to 



In this section we show in general that finite-step 
replica symmetry-breaking solutions for the fluid phase 
of the vortex lattice with a Gaussian disorder correlation 
function (10) are not stable. This will be shown irre- 
spective of the number of steps. We have shown this in 
the case of one-step replica symmetry breaking in Sec- 
tion VI. From (45) and (58) we obtain in the case of a 
i?-step replica symmetry-breaking solution in the fluid 
phase 



A/v, 



knT 



E 



1 



1 



nil 



with 



4arcsinh 



.1/2 



(l-hx/4)i/2 



(CI) 



(C2) 



^|E^[5(A™,)-5(A™,,J] 



-.9(Afl,) 



(C3) 



Y,— 5(A„,^)-5(A„,^_J 
for ^ = 1, . . . , i? with 



= 2?(2B[A„J)-Zz 

(C5) 



^^--EE^f--^) 

X [5(A^J - 5(A^^_ J] > . (C6) 



One should compare this equation with (60) in the case 
of a one-step replica symmetry-breaking solution. Since 
Zi = we can use (C4) and (C5) similarly as in the 
derivation of (81) to obtain A(l; 2, 2) < irrespective of 
the number R of hierarchical steps for Ai > 0. Further- 
more, we obtain A(l;2,2) = at Ai = as was also 
the case in the one-step hierarchical symmetry-breaking 
case. 
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